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ABSTRACT 

Combining redshift and galaxy shape information offers new exciting ways of ex- 
ploiting the gravitational lensing effect for studying the large scales of the cosmos. 
One application is the three-dimensional reconstruction of the matter density distri- 
bution which is explored in this paper. We give a generalisation of an already known 
minimum- variance estimator of the 3-D matter density distribution that facilitates the 
combination of thin redshift slices of sources with samples of broad redshift distribu- 
tions for an optimal reconstruction; sources can be given individual statistical weights. 
We show how, in principle, intrinsic alignments of source ellipticities or shear/intrinsic 
alignment correlations can be accommodated, albeit these effects are not the focus of 
this paper. We describe an efficient and fast way to implement the estimator on a 
contemporary desktop computer. Analytic estimates for the noise and biases in the 
reconstruction are given. Some regularisation (Wiener filtering) of the estimator, ad- 
justable by a tuning parameter, is necessary to increase the signal-to-noise to a sensible 
level and to suppress oscillations in radial direction. This, however, introduces as side 
effect a systematic shift and stretch of structures in radial direction. This bias can be 
expressed in terms of a radial PSF, comprising the limitations of the reconstruction 
due to source shot-noise and the unavoidably broad lensing kernel. We conclude that 
a 3-D mass-density reconstruction on galaxy cluster scales (~ 1 Mpc) is feasible but, 
for foreseeable surveys, a map with a S/N > 3 threshold is limited to structures with 
M 2 oo >lx 1O 14 M /i _:l , or 7 x f 14 Mq^ 1 , at low to moderate redshifts (z = O.f 
or 0.6). However, we find that a heavily smoothed full-sky map of the very large- 
scale density field may also be possible as the S/N of reconstructed modes increases 
towards larger scales. Future improvements of the method may be obtained by includ- 
ing higher-order lensing information (flexion) which could also be implemented into 
our algorithm. 
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1 INTRODUCTION 



Since i ts detection 
200d; iKaiser et af 



about 
l200d ; 



a decade ago 
IVan W acrbck e et al 



(iBacon et al.l 



200d ; 



Wittman et al.M2000f) weak gravitational lensing, or cosmic 
shear - the small, coherent distortion of distant galaxy 
images due to the intervening large-scale matter distribu- 
tion in the cosmos - has become a well-established tool 
for studying the dark Universe (for recent reviews see : 
Hoekstra fc Jainll2008l; ISchneidedl20o5 iMunshi et al"1l2006l ; 
Van Waerbeke fc Mellied l2003h . The phenomena of gravi- 
tationally scattering of light by massive structures is well 
described by General Relativity. Hence lensing provides us 
with a direct probe of the mass-distribution in the Universe. 



The spectrum of applications for exploiting this pool of 
information includes the interpretation of correlations of the 
cosmic shear signal to constrain cosmological parameters, 
the cross-correlation of cosmic shear with galaxy positions 
for studying the typical matter environment of galaxies or 
galaxy biasing, and the fascinating possibility of mapping 
the large-scale mass distribution in the Universe. 

Provided that General Relativity itself is accurate 
( Hute rer fc Lindedl2007l :lu zan fc Bernardea"u1l200ll ) the sta- 
tistical analysis of the weak shear signal for cosmological 
parameter estimation provides us with ample indirect evi- 
dence that most of the matter of the Universe is not made 
up of familiar baryonic matter, but of so-called Dark Mat- 
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ter. Moreover the largest piece in the energy budget in the 
present Universe is occupied by the so-called Dark Energy 
which seems to be accelerating the expansion of the Uni- 
verse (see e.g. iPeacock et ai1l2006l ; lAlbrecht et"al]|2006l for 
a recent review). As the area, A, of weak lensing surveys 
grows, the statistical uncertainty on the parameters of this 
cosmological model will decrease as l/\f~A, allowing even 
General Relativity itself to be tested. The main issue for 
lensing surveys for extracting cosmological p arameters will 
be th e control of systematics at the 1%-level |Huterer et al.l 

l2ooeh . 

As well as statistical analysis, weak lensing can also 
be used to directly visualise t he m atter distribution in 
the Universe. iKaiser fc Squires! (|l993h first showed how to 
ma p the projected 2-D mass distribution fro m weak lens- 
ing. iTavlod |200ll '). lBacon fe Tavlorl (|2003h and lTavlor et al.1 
(2004) showed how to extend this to 3-D mapping by com- 
bining shear data with redshift distances (3-D weak lensing). 
Visualisation of the Dark Matter distribution has the clear 
advantage that one can directly compar e the matter distr i- 
bution with that of galaxies and gas (e.g. lClowe et al.l l2006). 

For a given resolution scale, mass mapping does not 
benefit f rom a larger area. The noise of a mass map s cales as 
<r € /y/N (jKaiser fc Sauireslll993l : iHu fc Keetonl I2002T ) where 
TV = n8 s is the surface number of lensing relevant galax- 
ies (sources) within a resolved area, 8 3 , the mean surface 
density of galaxies is n, and a E is the intrinsic random el- 
lipticity distribution of sources. It is unlikely that n and cr E 
will change compared to current surveys as dramatically as 
the survey area, while the uncertainty in redshifts must also 
be added for 3-D mass mapping. 

Ambitious survey campaigns are either under way 
(e.g. CFHTLS, RCS2, COSMOS), about to commence (e.g. 
Pan-STARRS, VST-KIDS, DES) to cover larger areas of sky, 
while surveys being planned for the near future (e.g. LSST, 
Euclid) will map the full sky to great depths and high res- 
olution. This will provide the lensing community with large 
galaxy catalogues with billions of galaxies including multi- 
colour information for photometric redshifts and accurate 
galaxy image measurements. These surveys mark the begin- 
ning of a new era for lensing as the availability of a large set 
of sources with redshifts and multi-colour data will allow 
us to accurately probe the evolution of mass-clustering as a 
function of cosmic time. 

This paper addresses in detail the question of how to 
optimally reconstruct the 3-D matter density solely based 
on a lensing catalogue with redshift information. Apply- 
ing Wiener filtering to cosmic shear follows in the tradi- 
tion of other successful cosmological appl ications, initially 



outlin ed for astrophysical applications by iRvbicki fc Press 



1992) and applied to galaxy redshift surveys ( Lahav et al 



ll994l;IZaroubi et al.lll995l; IWebster et al.lll997h . CMB maps 
( Bunn et al.1 Il994l; iBouchet et al. I ll999T ) and peculiar ve- 
locity fields (|Zaroubi et al.1 1 1999). Here we describe a lin- 
ear minimum-variance reconstruction operator which is a 
generalisation of th e Wiener filter method proposed by 
IHu fc Keetonl l|2002l '). iBacon fc Tavlorl (|2003T ) also proposed 
a Wiener filtering method, which was first applied in 



iTavlor et all (l2004h for a 3-D reconstru ction of the COMBO - 
17 survey (|Wolf et al.1 l2004h and in iMassev et ail (|2007f l 
to the COSMOS survey. However this was to reconstruct 
the gravitational potential, rather than the mass-density 



field, which improved the noise pr operties of the 3-D re - 
construction on clus ter scales (see IBacon fc Tavlorl [20031 ). 
IHu fc Keetonl (|2002T ) pointed out that on larger scales 
the si gnal-to-noise improv es all owing density-reconstru ction. 
Both lHu fc Keetonl (|2002T l and lBacon fc Tavlorl 120031 ) pro- 
posed reconstructing the 3-D convergence field in a first step, 
and then applying the Wiener filter along each individual 
line-of-sight, ignoring the correlation of errors between dif- 
ferent line-of-sights (radial mapmaking). In the application 
to the COMBO-17 and COSMOS surveys, all of the redshift 
bins where assumed to have the same noise level, which is 
sub-optimal. 

In this paper, the transition from 3-D shear to lens- 
plane density contrast is done in one step, along all line-of- 
sights simultaneously, including all correlations between dif- 
ferent lines-of-sight and lens planes. Another extension con- 
sists in introducing an additional parameter that regulates 
the impact of the Wiener regularisation. With the approach 
outlined here we show how all available redshift information, 
i.e. redshift slices and sub-samples with broad redshift distri- 
butions and individual statistical weights of sources, can be 
combined in one framework. The primary focus of this study 
will be the signal-to-noise and biases that are introduced in 
the 3-D reconstruction if we regularise the process with a 
prior. Regularisation will be unavoidable due to the poor 
signal-to-noise of a 3-D density map based on lensing alone. 
Besides the theoretical considerations, we also describe how 
the reconstruction algorithm can be efficiently implemented 
in practise. 

This paper is laid out as follows. We start off in Sect. 
[5] with a very general discussion of inverse estimators for 
linear problems. Sect. [3] continues by specifying the prob- 
lem of 3-D reconstructions in two steps: the convergence to- 
mography and the lens-plane matter density contrast from 
the cosmic shear tomography. Both problems fall into the 
category of linear inversion problems. Sect. [4] discusses the 
signal-to-noise of mass density reconstructions, the main fo- 
cus of this paper, and biases therein that are introduced by a 
prior. As prior we consider expected density correlations in 
radial or transverse direction. The practical implementation 
of the estimators, Appendix [S] requires some optimisation, 
because a naive implementation would force us to handle ex- 
tremely large matrices even on relatively small grids. Finally 
in Sect. [5] we apply the algorithm to a Active, simulated sur- 
vey with specification we should expect for a weak lensing 
surveys within the next twenty years. 

We use as fiducial cosmology a ACDM model (adiabatic 
fluctuations) with fi m = 0.24, tt b = 0.0416, Oa = 1 - !Jm 
and H = h 100 km s" 1 Mpc" 1 with h = 0.732. The normal- 
isation of the matter fluctuations within a sphere of radius 
8/i _1 Mpc at redshift zero is assumed to be as = 0.76. For 
the spectral index of the primordial matter power spectrum 
we use n s = 0.96. Th ese values are consis tent with the third- 
year WMAP results (|Spergel et al.ll2007f ). 



2 GENERAL LINEAR INVERSIONS 

2.1 Estimators for linear problems 

All problems discussed in this paper are linear inversion 
problems. For such problems, a vector of observables, d, is 
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known to be a linear transformation of a certain vector of 
underlying parameters, s: 

d=Rs + n, (1) 

where R is a linear projection from s-space to d-space. The 
vector n denotes a set of random variables that introduces 
noise to the observables. The statistical uncertainties of the 
observables are thought to be vanishing on average, (n) = 0, 
and may be a function of s. All involved matrices and vectors 
may be complex numbers, the dagger denotes the complex 
conjugate transpose. 

The inversion of Eq. [1] with respect to s for a given d 
can, in general, only be done in a probabilistic way because 
of a) the noise in the observation and b) because R is usu- 
ally not regular or not even a square matrix. One approach, 
the maximum-likelihood (ML) approach, uses the posterior 
likelihood function for s given the observation d and known 
probability distributions of the involved random variables 
to fin d the most likely solution for s (e.g. iKitaura fc EnBlir] 
l200Sh : 

P(s\d) = oc P(d\s)P(s) . (2) 

The Bayesian evidence P(d) is a constant with no interest 
for the scope of this paper, P(s) is a prior on s and P(d\s) 
the likelihood function of d given some parameters s. 

In a weaker approach to the inversion problem, one can 
demand to find a linear operator, Hd say, acting on the ob- 
servation d that minimises the variance of residuals between 
the guessed signal and true signal, {(s — Hd)(s — HdV), 
witho ut the strong assumption of P(s) and P(s\d) (Seliak 
1998). All that is now required, are the statistical means of 
s and n, both assumed to be zero in our case (or can be 
made zero by definition), and their covariances 
and S = {ss^}. Minimising the covariance of the residuals 
for all conceivable noise and signal configurations results in 
a minimum variance (MV) estimator: 

s M v = Hd = |^S _1 + R t N _1 Rj _1 R 1 ^ x d (3) 

= ^1 + SR t N _1 Rj _1 SR*N _1 d , (4) 

where 1 is a unity matrix of the dimensions of S. We call 
this filter H a Wiener filter. Here, we have explicitly assumed 
no correlation between signal and noise, i.e. (sn^) = (see 
Appendix [B}. 

If we would not like or cannot make strong assumptions 
on the signal covariance, we can assume it to have infinite 
variance, S _1 = 0, avoiding any preference in parameter 
space and only using the noise covariance as weight (inverse 
variance filter): 

siv = ^IM^r] 1 R f N x d . (5) 

In the Bayesian interpretation of the MV-estimator we 
would not have any prior on s. 

For multivariate Gaussian distributions P(s) and 
P(d\s) the maximum-likelihood and minimum-variance so- 
lution are identical, so that in the Gaussian case the most 
likely solution is the minimum-variance solution. For this 
paper, the reconstructed fields will be cosmological density 
fields smoothed to large scales which makes them roughly 
Gaussian fields. 



Commonly, the sizes of the vectors and matrices can 
become large so that multiplications between vectors and 
matrices , a Q(N 2 ) proce ss, or matrix inversion, a 0(N 3 ) 
process dPress et al.ll 19921 ). can take a very long time or re- 
quire much computer memory if done without optimisations. 
The last step in Eq. U ensures that we have one matrix in- 
version less. 

It may be that in the estimators Eq.|4]or, in particular, 
Eq.[S]the matrices in the square brackets are not regular and 
hence cannot be inverted. This happens in reconstructions 
of the lensing convergence due to the mass-sheet degeneracy. 
This problem can be solved by enforcing another regulari- 
sation that asserts the mean over all elements of smv to be 
exactly a particular constant so. This makes the minimisa- 
tion problem of finding H a minimisation w ith additional 
const raint s^s = so breaking the degeneracy |Hu fc Keetonl 
l2002j y 

Alternatively, we can assert s to be close to zero (or 
another constant value) confined by a Gaussian probability 
distribution as prior oc e _ 2 Asts ; A is a constant constrain- 
ing the allowed regime of s^s. This is achieved by adding 
another equation to the minimisation of the aforementioned 
residual matri x for H and re sults in the so-called Tikhonov 
regularisation (|Wrightlll996l ) also breaking the degeneracy: 

s M v= [Al + RtN-'R] -1 RtN^'d. (6) 

Note that using the Wiener filter, and therefore a constraint 
on the allowed variance of s's via S, usually already does 
the job and breaks, in particular, the mass-sheet degeneracy. 

A fusion of a Wiener filter and a filter without prior, i.e. 
S _1 = 0, is introduced by a tuning parameter ( "a-tuning" ) , 
a £ [0, oo [, for the Wiener filter that changes the impact of 
the Wiener filter by scaling up the expected signal-to-noise 
in the data by S -1 t— > aS _1 (Saskatoon filter: lTegmark1ll997l : 
iTegmark etalll 19971 ): 

s M v= [al + SRtlNT'R]" 1 SR^^d . (7) 

For a = 1, the Saskatoon filter is identical to the normal 
Wiener filter, and for a — we would apply the filter with- 
out any prior. The Saskatoon filter will be of central im- 
portance for the reconstruction technique outlined in this 
paper. 

2.2 Covariance 

The (noise) covariance of the estimator Eq. [7] is found by 
working out 

Cov(smv) = (smvSmv) _ < s mv)(s[ iv ) , (8) 

with (...) being the ensemble average over all realisations 
of the noise n but leaving the embedded signal unchanged. 
This results in the expression 

Cov(smv) = WXW f . (9) 

Here, X= [R^N : R] 1 is the covariance of smv for the 
no-prior case and W = S[S + aX] _1 the Wiener filter in s- 
space. The covariance matrix has two limiting regimes. In 
the regime where the data is noise dominated, SX _1 <g; 1, it 
is approximately SX _1 S. On the other hand, for data where 
the noise is small compared to the signal, |SX _1 | >• 1, one 
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finds X, which is exactly the covariance of the MV-estimator 
without any Wiener regularisation. 

The power of the signal that gets through the filter (on 

average) , 



Pow(smv) = (smvSmv)' - < S MV)'(S^ V }' , 



(10) 



obtained by taking the ensemble average (...)' over all pos- 
sible signal realisations while setting the noise n = 0, i^j 



Pow(smv) = WSW f 



(11) 



In the signal dominated regime, this is roughly S, but in the 
noise dominated regime it is well described by SX _1 SX _1 S. 
This shows that a Wiener filter suppresses the signal if the 
signal-to-noise becomes low. 



2.3 Bias 

We briefly address here the question if the MV-estimators 
are unbiased estimators, i.e. whether (smv)' = s. Using the 
above equations and (d)' = Rs we find 

(smv)' = Ws (12) 

for the Wiener filter. For sufficiently small noise, N -1 ^ oo 
the Wiener filer asymptotically approaches XX" 1 , which is 
the identity matrix if R^N X R = X -1 is invertible and un- 
defined otherwise. We may call problems in which X -1 is 
invertible well-posed. Note that in the forgoing discussion 
of the asymptotic estimator covariances we assumed a well- 
posed problem. 

Therefore, for well-posed problems the Wiener filter is 
asymptotically unbiased for an increasingly higher signal- 
to-noise, and the non- Wiener filter is always unbiased (no 
prior). As we will see later, the Wiener filter applies a 
smoothing to the reconstruction by rescaling of low S/N- 
modes to some degree which explains why it is not an unbi- 
ased estimator in general. 

Non- well posed problems loose information when apply- 
ing R to the original signal. Therefore, even with no noise one 
cannot recover the exact signal. Those problems can be made 
treatable, by introducing the Tikhonov regularisation, if we 
make X -1 + Al invertible. This, however, will also to some 
extent bias the estimator. One finds W = (X -1 + A1) _1 X _1 
is the bias matrix in this case. 

For a Wiener filter an original signal in one of the com- 
ponents of s is, on average, in the reconstruction rescaled 
and spread over the other vector components smv- The ith 
column of the (bias) matrix W tells us how the original sig- 
nal Si is spread over the other vector elements of smv in the 
reconstruction. 

The apparently obvious thing to remove the bias in the 
estimate would be to multiply the MV-estimator by the in- 
verse, if it exists, of the bias matrix: 



W 



smv 



= XR f N _1 (Rs + n) = s + XR f N _1 n 



(13) 



This, however, yields exactly the estimator Eq. [5] removing 
the benefits of the Wiener filter. The noise covariance is now 
X again as can be seen from the second term in the r.h.s. of 
the last equation. 

1 It is assumed that the covariance of the true signal is indeed 
identical to the Wiener prior, i.e. (ss^)' = S. 



The Saskatoon filter, Eq. [7] is a compromise between 
the advantages of the Wiener filter, Eq. [4] (improved signal- 
to- noise but biased), and the no-prior filter, Eq. (noisy but 
unbiased). However, less Wiener filtering means less smooth- 
ing in the, then noisier, reconstruction which here will re- 
quir e an add i tional (non-optimal) smoothing after filtering. 

IZaroubil (|2002| ) derive another linear filter, a specific 
quadratic estimator, for unbiased reconstruction of the 
large-scale matter distribution from galaxy redshift surveys. 
Principally, this filter cold be employed in the context of this 
p aper as w e ll. Ho wever, as it turns out, the UMV-estimator 
of IZaroubil (|2002T ) is identical to Eq. [7] with setting N = 1 
(equal weighing of all data points) and a = 0, which is 
the no-prior filter, Eq. [5] with N = 1 and is therefore auto- 
matically included in the following discussion. Moreover, as 
already pointed out, the no- prior filter is singular for lensing 
applications so that some regularisation, such as a > or 
A > 0, would be required for the UMV-filter, then being a 
biased as well. 



3 RECONSTRUCTION METHOD 

3.1 Lensing convergence tomography 

Let us see how this fits into the context of three-dimensional 
mass reconstructions based on c osmic shear tomography. 
In weak gravitational lensing (e.g. iBartelmann fc Schneider] 
l200l]). the observables are complex ellipticities, E^ 1 with k 
being a galaxy index, of galaxy sources that, for our pur- 
poses, belong to one out of i = 1 . . . N z different redshift 
bins. The source ellipticity is, in the weak lensing regime, 
an unbiased estimator of the cosmic shear, 7, at the posi- 
tion, Ok, of the source. An angular position on the sky is 
denoted by 9 which, for convenience, is a two-dimensional 
Cartesian vector expressed as complex number. 
For the following, we assume that the vector 



7 



(0 



.(0 Ji) 



r 



(14) 



contains the complex ellipticities binned on a grid, N s is the 
number of grid points. Only sources belonging to the same 
redshift bin are binned together onto an angular grid. All 
source ellipticities, E^\ contained inside the fcth grid pixel 
and the same redshift bin, i, are combined to one average 
ellipticity 



(i) E (i) 



E 



(0 



(15) 



where E- and wn are the ellipticity and statistical weight 
of the sources associated with the fcth grid pixel, respec- 
tively. If there are no sources inside the pixel or the sum of 



all w 



(0 



0. We will further assume that 



, is vanishing, set e k 
sources of all redshift bins are binned onto the same angular 
grid so that same elements el and eP of different bins i 
and j lie along the same line-of-sight, Ok- By A we denote 
the solid angle of the grid pixels. Naturally, the number of 
grid points is identical for all redshift bins. 

The affiliation of a source to the ith redshift bin means 
that we have prior knowledge on its redshift. We denote this 
prior probability on the redshift, z, by p^\z). Although we 
are using the term "redshift bin" for pt\z) in this context, 
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we would like to point out that the "bins" may be arbitrary, 
including overlapping, distributions of sources in redshift. 
It may be conceivable, for example, that only a small frac- 
tion of sources have redshift information and can thereby 
be ordered into non-overlapping z-bins. The cosmic shear 
information of all the other sources could then be added as 
well by putting them into one "bin" with a wide redshift 
distribution, which is overlapping with all the other bins. 
This way, we can combine all available shear catalogues to 
acquire the best possible reconstruction^ 

In the weak lensing regime, the lensing convergence 
- also arranged as JV g - dimensional vector of convergences at 
position Ok analogues to 7^' - is related to the source ellip- 
ticities by a linear transform ation P-, K , a complex 7V g x iV g 
matrix (|Hu fc KeetorJl2002f ): 



7 



(0 



(16) 



where n\ is a vector of random intrinsic s ource ellipticities. 
The l inear transformation P 7K is for k ^ I IjKaiser fc Squired 
1 19931 ) 



A 1 



~tK,\kl 



(17) 



where 0% t is the complex conjugate of the separation vector 
61 - 6 k . For fc = I, we set [P 1K ] k i = 0. 

Note that, so far, we have carefully separated the dif- 
ferent redshift bins from each other by using the superscript 
"M» p or eacn bin i there is one Eq. 1161 We can combine 
everything to one compact system of equations, 



7 = P-ynK + n-) 



(18) 



by putting together all information into the vectors: 

7 = (7 (1) ,7 (2) ,...,7 (JV ^) , (19) 



AN,) 



)■ 



AN,) 



= diag{P 7fJ ,P 7 



(20) 

(21) 
(22) 



The round brackets, grouping together the vector argu- 
ments, should be understood as big vectors obtained by pil- 
ing up all embraced vectors on top of each other. Employing 
the linear relation between shear and convergence and plug- 
ging them into the Saskatoon filter, Eq. [7] with R = P 7 «, 
S = (kk')' and N = (n 7 n 7 ) gives an estimator, kmv, for 
the lensing convergence seen by the sources within the var- 
ious redshift bins. 

This estimator is ideally, only for exactly vanishing B- 
modes in the cosmic shear signal, a purely real vector. The 
imaginary part of kmv is at best pure noise, but contains in 
general a significant signal owing to a possible B-mode signal 
in the cosmic shear field. Thus, the imaginary part of kmv 
can be used as parallel B-mode reconstruction, whereas the 
real part is the sought E-mode of the convergence tomog- 
raphy. E-mode and B-mode are interchangeable by rotating 
the source ellipticities by 45°. 

2 Note that the following assumes that the noise covariance is 
diagonal. If we have the same galaxies belonging to different bins 
this is only a good approximation if their fraction is small. If not, 
we have to account for this in the noise covariance. 



For the statistical uncertainty, needed for N, of the 
source ellipticity attached to the fcth grid pixel and source 
bin i, we take 



k 



with 



H'T 



(i) 



,(0 



.(0 



1 

£ HI 



(23) 

(24) 
(25) 



Here the sum for the overall ellipticity variance, [a' 1 '] 2 , runs 
over all source ellipticities belonging to the ith source bin, 
and the sum for the effective number of sources inside the 
fcth pixel, only over the statistical weights of sources 

belonging to that pixel. For pixels containing no sources, 
we set CTe = (Too with CToo = 10 3 . The noise covariance is a 
diagonal matrix, 8^ and denote Kronecker symbols, since 
cross-correlations of the intrinsic ellipticities of sources are 
neglected. 

Although it is assumed for the scope of this paper that 
intrinsic ellipticities of sources are uncorrelated to any other 
relevant quantity, it may be possible to accommodate also 
cases where this is not true. For example, if we expect in- 
trinsic alignments between sources of the sam e or different z- 
bins, i.e. correlations of intrinsic ellipticities (|Heavens et al.l 
|2000| ; ICroft fc Metzled |2000| ; ICrittenden et al.l l200ll ). This 
case can be treated by adding the two-point correlation of 
intrinsic alignments, (8) as additional noise contaminant 
to N: 



rK jK 
; °kl°ij ' 



,eoi2 



(26) 



Another possible known contaminant of lensing ap- 
plications is the possible correlation betw een the intrinsic 
align ment of sources and the lensing signal dHirata fc Seliakl 
[2004J). This would give rise to a correlation between noise 
and signal, (sn^) 7^ 0, which was excluded here from the 
very beginning of the derivation of the linear filters and is 
being assumed to be negligible for this work. We defer the 
treatment of this contaminant to a future paper but already 
outline in Appendix [B] how this can in principle be included 
in the estimator. 

The signal covariance, S, is the two-point correlation 
function of the lensing convergence at the separation of two 
grid points, which is 



■|6 (i) ]*[K 0) ]j 



)' = (7 (i) (e k )[7 0) W)' = e? i) (f«). (27) 

where £ ^ is a two-point correlation funct ion of the cosmic 
shear (e.g. iBartelmann fc Schneider! 120011 ) between the ith 
and jr'th source bin, 7^(0fc) and ^\9i). 



3.2 Density contrast on lens planes 

We can go one step further than the lensing convergence 
by recognising that the convergence k' 1 ' itself, as seen in 
the sources of the zth redshift bin, is within the framework 
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ens plane redshift z 

Figure 1. Lensing efficiency, Qij, for j = 1 . . . 15 lens- 
planes, equally spaced between redshifts z = 0, 1.5(Az = 0.1), 
and i = 1 ... 30 source redshift bins, equally spaced between 
z = 0, 1.5(Az = 0.05). The numbers labelling some of the lines de- 
note the mean source redshift; the source redshift increases from 
left to right. The galaxy distribution within the redshift bins is 
assumed to be flat. 



of general relativity a linear transformation of the matter 
density contrast, S, along the line-of-sight: 



W (l) (w) 



3#o 

2c 2 



(28) 



dw 



W (i) (w)f K {w) 



dw 



a(w) 
, f K (w' - w) 
/ K (0 



5 (f K (w)9k,w) 
dz~ 



dw 



(29) 



= z(w') 



The function /k (w) is the angular diameter distance as func- 
tion of the comoving radial distance w and the curvature, 
K, of the fiducial cosmological model: 

f K- 1 ^ 2 sm(K^ 2 w) 
/k(w) = < w 

{ {-K)- 1/2 smh[(-K) 1/2 w] 



(K > 0) 
(K = 0) 
(K < 0) 



(30) 



Furthermore, we denote the vacuum speed of light and the 
cosmological scale factor by c and a(w), respectively. For the 
density contrast, we have chosen the comoving coordinate 
frame such that the first argument x± in 5(x±,xn) denotes 
a position in a plane perpendicular to the line-of-sight, and 
the second argument, x\\ , the distance of the plane along the 
line-of-sight. Note that we are assuming a flat sky. 

For a practical approximation of Eq. 1281 we postulate 5 
to be constant over I — 1 . . . N\ p comoving distance ranges 
x\\ e [wi,w l+1 ]: 



[A 



3Hj 

2c 2 



w { ) (w /k w 
dw -— 



(31) 



(32) 



This means, we assume that the 3-D density contrast stays 
constant along the line-of-sight as long as we stay inside 
the same matter slice; 6^ is here the average 3-D density 
contrast, and not the total, over the width of the lens-slice 



integrated 3-D contrast as can also be found in the literature. 
The coefficients Qu (lensing efficiency) define the response 
of the convergence tomography to the matter density on 
the individual lens-planes. A mean density contrast 5, over 
the radial width of the ith slice, induces on average shear 
polarisation of sources belonging to the jth redshift bin of 
amplitude QjiS. The coefficients Qy are independent of the 
line-of-sight direction, Fig. [1] displays an example for Qu. 

Due to the approximation in Eq. 1311 the mean matter 
density contrast within a slice is: 



r(0 



(w) 



Q+ 1 dw pM (w)S (/k (w)6 k , w) 

W {i \w)f K {w) 
a(w) 



(33) 



(34) 



hence a mean weighed with p^'(w). For relatively narrow 
lens-planes, the weights p^*' (w) will be roughly constant over 
the slice width so that they can overall be approximated by 
top-hat functions with width Awi = Wi+i —Wi. This is what 
we will assume for the remainder of the paper. 

As before, we can cast the relation between mean matter 
density and convergence tomography into a compact nota- 
tion by writing: 



k = Q<5 

where 
5 = 

Q = 



) (1 \S^ 2 \ 

( lQn 
IQ21 

\ 1Qn z i 



IQ12 • 



lQ2JV lp 



(35) 
(36) 
(37) 



and where 1 is the iV g x N g unity matrix. The relation be- 
tween the observable source ellipticities and the lens-plane 
density contrast is therefore: 



7 = P 7K Q# + n 1 



(38) 



As before for the lensing convergence tomography, an esti- 
mator for the average density contrast, <5mv, on the lens- 
planes is given by Eq. [7] with R = P 7re Q, N = (n 7 n 7 ) and 
S = Si = (SS 1 )' . For an efficient, practical implementation 
of the filter see Appendix IA1I 

Analogous to the estimator kmv, Smv will be a complex 
vector although the density contrast is by definition a real 
number. The real part of Smv will be the mass reconstruc- 
tion due to the cosmic shear E-mode, whereas the imaginary 
part has to be attributed to the cosmic shear B-mode. 



3.2.1 Transverse Wiener filter 

How can a prior on the signal covariance Sa be computed 
for a given fiducial cosmological model? By exploiting the 
Limb er equation in Fourier space (iBartclm ann fc Schneider] 
l200ll ). the relation between 3-D matter density and pro- 
jected density can be used to work out the second-order 



correlations of <5 



(0 X U) 



W. 



S 'S. 



Ji-K (i) 



(39) 
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6 in [ARCMIN] 

Figure 2. Expected projected two-point correlation, u)('\0), of 
the density contrast situated inside i = 1 . . . 9 equally spaced, 
Az = 0.1, lens-planes with redshift limits between 2 = 0... 0.9. 
The mean redshift of the lens-planes decreases from the bottom 
to the top line. The correlation functions have been smoothed for 
a grid with pixel scale 1'. 



1 f'°° 

± d££P^(£)\F(Q s l)\ 2 J (l6), (40) 
1 H+ 1 dw ( £ V 

(a^) 2 L [/kwi 2 ^Um'T ' 



P^(i) = 



Here P3d(k,w) is the fluctuation power spectrum of the 3- 
D matter distribution for a spatial scale k at a comoving 
distan ce w. For the scop e of this paper, the approxima- 
tion of ISmith et "all (|2003T ) will be used. The function P 4 W (I) 
quantifies the fluctuation power spectrum of (8), the pro- 
jected density contrast on the lens-planes. Note that cross- 
correlations between 5^''s belonging to different lens-planes 
vanish in the Limber approximation. Fig. Ogives an exam- 
ple for a/ l ^(#). The window kernel F(£) accounts for the 
binning on the grid, it is the Fourier transform of a grid 
pixel, which we assume as an approximation to be a circular 
top-hat function, A = 7r6> 2 , with radius Q B , hence: 

2Ji(e s ^) 



F(O s £) 



(42) 



J n {x) is the Bessel function of first kind. 

We call this prior on the correlations of density contrasts 
in different directions a transverse Wiener filter. 



3.2.2 Radial Wiener filter 

iHu fc Keetonl (I2002TI discuss an other Wiener filter, also ap- 
plied in iBacon fc Taylor! (|2003l ) for the 3-D potential, that 
puts priors only on the correlation between structures (av- 
erage density contrast inside pixel needle light cones) in the 
radial direction, i.e. 



where 
F(k) 



Ski 



(2tt) 3 



F(k)F*(k)P^\k) 



3 W) 



(fc) 



ik n wj sin(fc||Aw a /2) 

\J Ps(k, Wi)Ps(k, Wj) 
k«+k 2 ± . 



F(k ± w t O s ) 



(43) 



(44) 

(45) 
(46) 



The solid angle S denotes the radius of a circular pixel 
in radians and Wi = Wi + Awi/2 is the comoving (radial) 
distance of the centre of the ith lens-plane. 

In the total absence of correlations between pixels of 
different line-of-sights, the angular power spectra are flat 
power spectra (white spectra), which for a circular top-hat 
smoothing window equal^f] 



(47) 



We use this prior for the radial Wiener filter. Obviously, the 
amplitude of the flat signal power depends on the radius, O s , 
of the pixels since Eq. [43] is also smoothing density fluctu- 
ations in transverse direction. The following discussion will 
adopt <d s — 1'. This may appear too specific, but a change 
of the radius will only rescale the overall amplitude of Pij 
which, in turn, will solely act as a rescaling of the tuning 
parameter a. The essential feature of the radial filter, how- 
ever, namely its scale-independence, will remain unaffected 
by the choice of O s . 

We find that Py is close to diagonal with negative off- 
diagonal elements roughly two orders of magnitude or more 
smaller than the diagonal elements. Thus, the main differ- 
ence between our radial and transverse Wiener filter is the 
scale-independence of the radial filter, the off-diagonal el- 
ements are not of much importance which we checked by 
setting them to zero. 

The reader may be reminded that the two different 
Wiener priors discussed here - the radial and the transverse 
filter - differ only in the way assumptions about correlations 
between pixels on the lens planes are made. The terms do 
not mean that radial and transverse modes are reconstructed 
separately. All modes are determined simultaneously during 
the process of reconstruction, no decomposition into radial 
and transverse modes is done. 



3.3 Fourier space reconstructions 

For a discussion of the noise properties of the 3-D mass re- 
construction in the following section a Fourier representation 
of the problem is extremely helpful. Moreover, it allows with 
some idealistic assumptions to perform a quick-and-dirty re- 
construction in Fourier space which is of practical interest, 
see Appendix I A2I 

An idealised survey is considered. Such a survey has 
either periodic boundary conditions or is infinite in extent 
(flat sky), and the statistical noise ny is homogeneous for 
all source redshift bins. The latter means that the statistical 
uncertainty in the shear estimate is equal for each grid cell 
of the same redshift bin. 

In this case the matter density contrast, and all other 
fields, can equivalently be expressed in terms of its angular 
modes 



8(9) 



(27T) 



■8(1) e +1 



(48) 



3 As a related example, take the intrinsic shear ellipticity noise 
which has the fiat power a 2 n~ 1 . The intrinsic noise inside a cir- 
cular pixel is then (T 2 /fi7rG 2 = a 2 N~ 1 , where N is the mean 
number of sources inside a pixel. Hence, we would expect the in- 
trinsic noise variance to go down as <x 1/vJV (Poisson shot noise). 
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S(£) = / d 2 GS(G)e- 1 



i/2(0i*+e*£) 



(49) 



which defines the density contrast on a typical scale 2tv/\£\. 
The advantage is that now only angular modes of the same 
£ couple to produce the angular modes of the shear tomog- 
raphy because the conversion from lensing convergence to 
shear is a convolution, 



(£)=D(l)J2Qij~5 U) (*) + tiy\t) 



(50) 



3=1 



where D(£) = £/£* (|Kaiser fc Squireslll993l ) and £ = £i + i£ 2 
with l\ and £2 being the Cartesian components of the angu- 
lar mode. 

Devising the filter Eq.[7]of Sect.[2]one gets for the MV- 
density mode 



(51) 



Note that S and 7 (the Fourier transform of the grid- 
ded shear catalogues) are written here in a compact vec- 
torial form with mere N\ p /N z elements, grouping all lens- 
planes/source bins together, as outlined in the foregoing sec- 
tions; a is the usual tuning parameter. Here we furthermore 
use 



Ma 



(52) 



which - for the sake of a simplistic noise model used in 
the next section - assumes a shot-noise term onljQ; n t is 
the mean number density of sources belonging to the ith 
redshift bin. 

The transverse Wiener filter constraints the expected 
signal power by 

Al) (\£\) . 



(53) 



Noise and signal are subject to the same smoothing if we 
bin our data on a grid. Therefore, the effect of a smooth- 
ing window cancels out in the filter, since Sa and N appear 
always as ratios in the filter. 

The (expected) signal Ss is a diagonal matrix for trans- 
verse filtering in the Fourier space representation. This may 
change, if we choose to apply radial Wiener filtering, where 
only correlations between matter densities along the same 
line-of-sight are regularised. Since correlations between dif- 
ferent directions are neglected, the power spectra used for 
regularisation are flat power spectra (Sect. [3~2.2|l . indepen- 
dent of £: 



P 



(54) 



that is thought to be a linear, invertible transformation of 
S. 

Such cases could be (lens-plane-) smoothed density re- 
constructions, for which F' 1 ' would be a (invertible) smooth- 
ing operator acting on 5^ , or the 2-D gravitational potential 
on the lens-planes w hich can be pictured a s a special type 
of smoothing of d (i) (|Bacon fc Tavlorll2003T l: 



4>{0) = ~ I dV<5(0')m|0'-0| 
fin 



(56) 



or in a form for a discrete grid assuming that the grid sizes 
for <f> and S are equal: 



F (l) ] fci = - In \0i-0k 



(57) 



The density contrast is constant over the size of one (round) 
pixel. For 61 = 6 k one finds therefore [F* l '] fefc = — A/2-k. 
Note that Eq. [56] gives only one possible solution for the 
2-D gravitational potential because the potential for a given 
8^' (8) is not unique. 

The new, smoothed </>, is related to the data, d, by 



d 
F 



RF- 1 d) + n 



diag{F (1) ,F (2) ,...,F 



CN lP ) j 



(58) 
(59) 



Following the arguments in Sect. [2] it is easy to show that 
the MV-estimator for d> is 



F5a 



(60) 



This means, once we have got the MV-solution of S we can 
simply apply F to it in order to acquire the MV-solution 
of <j>. We do not have to go through the process of an- 
other full MV-reconstruction, if we want to further smooth a 
MV-solution. That is, provided the smoothing is invertible. 
Moreover, the covariance of 4>mv ls 



Cov(0n 



F t Cov( < 5i, 



(61) 



BIAS AND NOISE PROPERTIES OF 
RECONSTRUCTIONS 



As shown in iHu fc Keetonl (|2002l ) a mass reconstruction, 
even with as many sources as 100 arcmin -2 and with heavy 
smoothing, is very noisy. Therefore, an additional regular- 
isation of the reconstruction, by means of Wiener filtering 
for instance, is an absolute necessity. With regularisation, 
however, a careful analysis of the remaining statistical un- 
certainties and biases in the reconstruction will be required. 
This will be the focus of this section. 



3.4 Gravitational potential on lens planes 

We may ask the question how the MV-estimator of S is 
related to the MV-estimator of 

</> = FS (55) 

4 In the presence of correlations between intrinsic source ellip- 
ticities, the noise covariance needs additionally to be offset by 
[N]jj 1— > [N]y + Pfo(£), where P^\t) is the power spectrum of 
the intrinsic correlations between the ith and j'th source bin as 
function of scale \t\. 



4.1 Average signal-to-noise of reconstructions 

Devising the MV-filter in Fourier space, Eq. 1511 offers a con- 
venient way to estimate the covariance of the noise in the 
reconstruction. The following conclusions, however, will be 
based on the simplifying assumption that the noise pattern 
is homogeneous and that there are no gaps in the data, in 
particular the shear pattern of the full (infinite) flat sky is 
known. 

We start by considering the signal-to-noise for individ- 
ual ^-modes and assuming that the matter fluctuations of 
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the large-scale structure that provide the signal in the shear 
catalogue have cosmic average amplitudes. 

The covariance of the statistical uncertainties in Eq. [51] 
is according to Eq. [9] 



Nmv(q) = WXW 
with 

W = S s [aX + Si]" 1 and X = [C^N^C)]- 1 



(62) 



(63) 



The latter is the noise covariance without a Wiener prior or 
if the prior is unimportant. 

How much signal-to-noise can we expect on average for 
a 3-D mass reconstruction as function of angular scale? To 
answer this question, we compare, for a ACDM fiducial cos- 
mological model, the expected signal power as function of 
scale I, S, getting through the filter (Eq. [IT) , 



Smv(q) = WSW , 
to the noise power: 



|(a) of *» v (£) 



V 



WXW 



(64) 



(65) 



As expected density fluctuation power we take the prediction 
based upon the Limber equation. 

Since the noise covariance, N, scales as oc of/n, we can 
derive useful scaling relations for the signal and noise power 
with a fixed redshift distribution of sources. One finds that 
the tuning parameter inside the filter scales as 



(66) 



if we change the intrinsic shape variance from a t to a' t or 
the source number density from n to n'. Thus, increasing 
the number density to n' , for instance, yields a filter corre- 
sponding to n but with a different, smaller a. Similarly, we 
get: 



N M v(a) 
Smv(«) 



(a') 



Siuv(tt') 



so that the signal-to-noise of modes scales as 



(67) 
(68) 

(69) 



This property is useful for preparing plots of the estimated 
signal-to-noise spanning a wide range of fiducial surveys. 

Fig. [3] is a particular example of the signal-to-noise of 
reconstructed density modes using realistic survey parame- 
ters. For the source distribution with redshift a function 



p(z) 



(70) 



with zo = 0.57 was used; the mean redshift is z = 0.85. 
The distribution was truncated beyond a redshift of z — 1.5. 
The redshift distribution is somewhat shallow but not too 
unrealistic for contemporary surveys if we take into account 
that source redshifts are required. 

Obviously, even for a very optimistic source density 
of lOOarcmin -2 , Wiener filtering is utterly necessary to 
improve the signal-to-noise in the reconstructions. Simple 



transverse smoothing of the reconstruction obtained with- 
out Wiener prior, i.e. suppressing high-^ modes, does not 
suffice as even the lowest ^-modes have a signal-to-noise less 
than unity. 

The transverse or radial Wiener filter bring the signal 
up to a signal-to-noise of roughly unity, somewhat higher on 
very large angular scales. If we tune down the Wiener filter 
by lowering the a-parameter, we get a signal-to-noise that 
lies somewhere in between of the full Wiener filter and the 
very noisy reconstruction with no prior (not shown). These 
figures also demonstrate the attenuation of low-S /N modes 

by the Wiener filter, i.e. the ratio J [Smv {q)]u / 'P^\i). This 
tells us that the Wiener filter will, on average, not recover 
the original amplitude of a fluctuation mode. Instead, the 
amplitude will, depending on the redshift of the lens plane 
and the scale £, be attenuated. This is an unavoidable side- 
effect of Wiener filtering, which, on the other hand, thereby 
improves the signal-to-noise in the reconstruction. Another 
side-effect concerning the radial distribution of densities will 
be discussed in the next section. 

For the transverse filter especially, small scale modes are 
damped, which thereby acts as a low-pass filter performing 
an automatic smoothing on the lens-planes. In case of a ra- 
dial Wiener filter, the damping is roughly the same on all 
scales due to the scale-independence of the regularisation. 

Notice that in Fig. [3] the tuning parameters for the 
transverse, a, and radial filter, 0, are not completely equiva- 
lent, although they have the same effect: For our transverse 
filter, the theoretical signal power is identical to the expected 
signal power plugged into the filter as Si = S (both based on 
Limber's equation), whereas for the radial Wiener filter, the 
expected power S - still the same as before to have a com- 
parison assuming the same input signal - is different from 
Si used inside the Wiener filter (flat power spectrum). 

The improvement of the signal-to-noise due to Wiener 
filtering comes with another price that has to be paid. The 
statistical errors between the modes of different lens-planes, 
which are weakly correlated if no Wiener filtering is applied, 
become correlated. Fig. [3] shows the matrix 



[N M v(a)]i 



'NMv(a)]u[N M v(a)]j 



(71) 



for the previous survey parameters with 30 arcmin - source 
density. With no prior at all one finds that the errors be- 
tween different lens-planes quickly decorrelat e in a oscilla- 
tory m anner (bottom right). As pointed out in lrlu fc Keetonl 
(2002) this particular form of error correlations is owed to 
the finite difference approximation of a fourth-derivative in 
radial direction. 

By means of the Wiener filter, however, the errors be- 
come correlated over a wider range (the filter also smoothes 
in radial direction). Errors between different £'s remain un- 
corrected as they were. For a signal-to-noise close to one, 
the error correlations limit our ability to pin down the red- 
shift of structures along the line-of-sight, that is the redshift 
resolution of the reconstruction. The correlations have also 
a positive side: the radial oscillations become stretched out 
by the Wiener filter. 

The foregoing way of estimating the signal-to-noise in 
a reconstruction takes a quite pessimistic view. Namely, it 
presumes that the matter fluctuations of the structure we 
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Figure 3. Signal-to-noise ratio of reconstructed matter density modes as function of angular scale for three different lens-planes (solid 
line: z = 0.16, dashed: z = 0.76, dashed-dotted: z = 1.05) using transverse (left figure) or radial (right figure) Wiener filtering (top lines; 
j3 = 0.1 for the radial filter, a = 1.0 otherwise) and no prior (bottom lines; a = 0). The tuning parameters of the radial, /3, and the 
transverse filter, a, are not equivalent (see text). The reconstruction uses 15 equally spaced lens-planes, between 2 = 0... 1.5, and 30 
equally space source redshift bins, between z = 0. . . 1.5, with a source number density of 30arcmin -2 (thin lines: lOOarcmin - 2 ) and 

cr t = 0.3. The mean redshift of the sources is z = 0.85. The right sub-panels show the attenuation factor y [Smv ( a )\ a Z^s M ( = 1 f° r 
no damping) inflicted by the Wiener filter as function of scale. Plots for a = are identical for both the radial and the transverse filter 
and are therefore only shown in the left panel. 



1=2160 




Figure 4. Correlation of statistical errors of 5fjy(£) and <5^jy(£) 
for three distinct I (transverse Wiener filter) and no prior or a = 
(bottom right), which is independent of t. This plot assumes 
cosmic average matter fluctuations. The Wiener filter assumes 
a = 1 and a galaxy number density of 30arcmin -2 . 



would like to recover are actually of an amplitude we expect 
for a cosmic average, i.e. we expect |5(^)| 2 ~ Ps(t). In fact, in 
individual regions of space we will find matter concentration 
exceeding the average (large galaxy clusters) , or we will find 
less matter (voids). 

Let's quickly make an rough estimate of how much 



power above the average we may can expect. For a Gaussian 
random field - probably being a fair approximation of the 
mass density field on large scales and, thus, the regime we 
can realistically hope to reconstruct - the expected variance 
of the fluctuation power is (exploiting Wick's theorem): 



P(lf ) - ( P(£) ) = 2P(£) 



(72) 



Therefore, we may expect, owing to the natural variance in 
a Gaussian random field, to find regions where the signal 
power is actually up to four (2a) times the average fluctua- 
tion power, enhancing the signal-to-noise in Fig. [3] by a fac- 
tor of two - or even more in the non-linear regime. In Sect. 
14.41 we will work out the signal-to-noise of specific matter 
haloes. 



4.2 Response in Fourier space 

There will be a fundamental limit in recovering the radial 
matter distribution, even if we have peaks with a high signal- 
to- noise. As discussed in Sect.[2J a Wiener based reconstruc- 
tion is unavoidably biased, unless the signal-to-noise of the 
data is infinite. This bias is quantified by the bias matrix W 
that relates the original vector of £-mode coefficients to the 
(statistical average) vector of modes in the reconstruction. 
It is a convolution applied to the true matter distribution. 

Specifically, the ith column of W expresses how the 
Wiener reconstruction, on average, responds in the recon- 
structed density contrast field to a value <r J '(£) in the true 
matter distribution. Ideally, we would like to have a unity 
bias matrix (unbiased estimator), but we find that, in the 
reconstruction, the Fourier coefficient becomes rescaled (at- 
tenuation), shifted and spread out in the radial direction, 
contributing to other <5^Jy(£) with i ^ j. This response for a 
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Figure 5. Examples for the response of the Wiener filtered 
(transverse) matter reconstruction algorithm to peaks in angu- 
lar Fourier space of coherence scale £ = 800 in the true matter 
distribution at low, z ~ 0.05, and high, z ~ 1.05, redshift (down- 
ward arrows). The fiducial survey has the parameters as outlined 
in the caption of Fig. [T] assuming a source density of 30arcmin -2 
with a = 1 (thick lines). In the 3-D reconstruction, the localised 
peaks are shifted in redshift and smeared out, independent of the 
underlying signal, which limits our ability to reconstruct the mat- 
ter distribution in radial direction. The thin lines correspond to 
a Wiener filter with n = lOOarcmin -2 and a = 1. Increasing n 
or decreasing a reduces the shift and spread of the response. 



certain £-mode is demonstrated by Fig. [5] In this particular 
but realistic example, an angular mode of a lens-plane at 
redshift z = 1.05 (z = 0.05) is shifted toz~ 0.5 (z = 0.15) 
and smeared out by a z « 0.3 (la- width; a z fa 0.1) to adja- 
cent lens planes. The smearing of the peak is the reason for 
the strong correlation of modes of neighbouring lens-planes, 
which has been discussed in the forgoing section. 

Especially structures at the upper redshift end of the 
survey a prone to the shift and smearing. Both effects are 
solely a function of angular scale, £, larger coherent struc- 
tures (smaller £'s) are easier to resolve than smaller struc- 
tures. For small enough £'a all redshifts are eventually shifted 
to one redshift and similar response functions. 



4.3 Point spread function 

The ^-dependence of the response in the reconstruction as 
described above is hard to relate to the question how well we 
can recover the redshift matter clump as function of mass 
and redshift. Therefore, the following translates the Fourier 
picture of the reconstructions back to real (angular) space, 
where we further study the Wiener filter impact. 

Assume we would like to calculate the effect of the 
Wiener filter on a halo sitting at some redshift, or equiva- 
lently on some lens plane i, with projected density contrast 
profile 



= [ dwS sd (f K (w)6,w) 
Aw Jo 



(73) 



Furthermore, we consider for simplicity only profiles that 
are rotationally invariant, i.e. 



A single "point" would be just a Dirac delta function in 
direction of 0. 

The prefactor Aw denotes the range of comoving dis- 
tances spanned by the matter slice the halo is located in. In 
our framework, &2d(0) is the average 3-D density contrast, 
constant inside the slice represented by the lens-plane, and 
not the integrated 3-D density contrast over the width of a 
slice, see Sect. 13.21 

If we place the halo at the centre of our coordinate 
system, — 0, and smooth the overdensity over the solid 
angle, A = 7r€> 2 , of the circular pixel with radius O s (our 
grid binning) centred on 6 = 0, 

8(e s ) = ^Jj 2 &8 2d (&) , (75) 

we would find on average after Wiener filtering, performing 
the integration over all modes £, a smoothed over-density of 



6(& s ,a) 



(2tt) 2 J R 2 
1 

2^ 



d 2 £F(£Q s )\N6(£) 



d££F(£Q s )\N5(£) 



(76) 
(77) 



where 5(1) is the 2-D Fourier transform of the radially sym- 
metric halo profile: 



6(£) = 2tt / d6 6 J (£6) 6(6) 



(78) 



The window F(x) defines the smoothing window (also radi- 
ally symmetric) of our lens-plane pixels, which for circular 
pixels is denned in Eq. 1421 

For the scope of this analysis, the original halo signal, 
5(£), is a vector of mode amplitudes on the lens-planes, being 
zero except for the plane on which the halo is located. After 
filtering, l.h.s. of Eq. 1761 the signal will be spread along the 
line-of-sight, i.e. spread over the output vector o(O a ,a), as 
already seen for the Fourier modes, Fig. [5] Moreover, as 
everything is linear, the response of a sum of haloes is just 
the sum of the individual responses. 

There is also a transverse (spread inside same lens- 
plane) point spread function (PSF) associated with the fil- 
ter, which can be evaluated at separation A9 from the halo 
centre by offsetting the position of 6(0) in Fourier space, 



/2(ei*+e*i) 



6(1) i-> 6(1) e [ 

and reassessing Eq. [75] (for a radially symmetric profile) 



(79) 



6(Q B ,a,A6) = — 



d££F (£e B ) Jo (I AO) \NS(£) 



(80) 



5 2d (0) = 5(\0\) 



(74) 



Obviously, the transverse PSF is symmetric (only a function 
of the modulus of A.6) and will, therefore, not bias the an- 
gular position of a halo. Moreover, usually a strongly down- 
tuned Wiener filter will be used, a < 1, which will make 
the transverse spread relatively small. As a matter of fact, 
we will have to apply an additional transverse smoothing to 
get a reasonable signal-to-noise data cube. Conversely, the 
radial PSF will be of more concern. 

Fig. [5] gives a toy model example of the response of the 
Wiener filter to some halo with profile 5(6, Zh), sitting at 
a singular redshift Zh — 0.5. If the response is a direction- 
independent function C(z^,z) - in a toy model fashion a 
Gaussian with mean z^ and width <r z = 0.2 - then we will 
find in the reconstruction the 3-D map 5(6, Zh)C(zu, z). The 



onno D Ac l\/ri\Tr> A c nr\n [Thiol 



12 Simon, Taylor & Hartlap 




Figure 6. Toy model example of the response to a halo with some 
density profile sitting at a redshift of z = 0.5, originally with no 
radial extension. Plotted are four different iso-density surfaces 
with decreasing iso-densities, growing bigger in size for smaller 
iso-density values. Clearly, the surfaces become stretched - the 
more, the smaller the iso-density - in radial direction in a cigar- 
like manner (we loose redshift information) due to the response of 
the Wiener filter. As shown in Fig- [8] the effect becomes smaller 
for smaller a (less Wiener prior) which, however, will increase 
the noise in the reconstruction. The iso-surfaces responding to 
structures situated at the redshift boundaries of the survey will 
be radially stretched in a asymmetric manner (not shown). 



general effect is that the iso-density surfaces of the profile are 
stretched radially. Although all 8(6, Zh) are being stretched 
with the same C(z^,z), surfaces with smaller iso-densities 
are stretched more than higher iso-density values. Qualita- 
tively, this toy model describes quite well what we find in 
the real reconstructions. Note that the surfaces are centred 
about the true halo redshift since no z-shift due to the re- 
sponse C(zh,z) was assumed for this illustration. 

Under the idealised circumstances the radial point 
spread function will be the same in all directions. Since, 
in reality, we will have a varying source number density as 
function of direction and a finite survey area we can expect a 
more complex, direction dependent response of the Wiener 
filter which accounts for inhomogeneous source ellipticity 
noise patterns. 

Assuming as a lowest-order approximation a homoge- 
neous noise pattern of the sources, we can work out the 
expected radial PSF of the Wiener filter, i.e. the radial con- 
volution C, by means of Eq. 1771 

i foo . 

C = 7T d£e\F{£Q B )\ 2 \N . (81) 

Here, an original pixel on the lens plane has the same shape 
and size as a reconstructed pixel, i.e. 5(£) oc F(£Q S ). A pixel 
on the ith lens plane is spread over the other lens planes 
according to Ce^, is a unity vector non-zero only at the 
original position of the pixel, or put in another way: dj is 
the signal-leakage of pixels on the ith lens-plane into pixels 
on the jth lens-plane. In an ideal, unbiased reconstruction, 
W = 1, the matrix kernel of the integral would be the unity 
matrix so that 

C = — / AU\F(£Q S )\ Z \ = 1 (82) 
2?r Jo 



for a properly normalised pixel window function F. Thus, 
the pixelised reconstruction would yield exactly the original 
matter distribution on the lens-planes smoothed with the 
pixel window. In reality, however, we will find a situation 
more like the one shown in Fig. [7] 

Now, Fig. [5]uses the results of the previous discussion to 
quantify the quality of a 3-D reconstruction for a particular- 
fiducial cosmology and a realistic fiducial survey employing 
a Wiener filtering scheme with a-tuning. The plots can be 
scaled to different cr e or n as explained in the caption. We 
can see in the plots that, the shift and spread of the signal 
in the reconstruction due to Wiener filtering decreases if the 
filter is tuned down, which, on the other hand, will decrease 
the signal-to-noise. We also find that the redshifts of objects 
near the "edges" of the survey (here: z ~ and z ~ 1.0) 
are mostly biased, high redshift objects appear to be more 
affected. The bias tends to shift low-z objects to apparently 
higher redshifts, while high-z objects are moved in direction 
of small redshifts. 

In analogy to these plots, Fig.|5]shows the same param- 
eters but this time obtained with a radial Wiener filter. We 
find a similar behaviour compared to the transverse filter, 
keeping in mind that the tuning parameters of the radial and 
transverse filter are not comparable: an (3 of the radial filter 
ten times smaller than the a of the transverse filter yields 
roughly the same signal-to-noise (see next section). The ex- 
pected signal based on the radial filter gives, in contrast to 
the transverse our filter, a different prediction as the "true" 
signal which is calculated from Limber's equation. Looking 
more closely, one finds that a radial filter might be used to 
slightly decrease the 2-shift due to the PSF ((3 ~ 10 -3 com- 
pared to a ~ 10~ 2 ), although the spread in radial direction 
may be larger. As this spread is partly into the negative 
regime, see right Fig. [7] it may not become immediately ob- 
vious but can have the effect of diminishing structures along 
the line-of-sight - or generating structures as a response to 
an underdensity. 

We noticed during the course of this work that a very 
simple radial prior that gives all lens-planes the same prior, 
Pij = 7T0g, gives often somewhat less biased signal-to-noise 
maps so that it is worthwhile to compare in a concrete case 
the benchmarks of this simple radial filter to the previous 
radial filter. 

As already pointed out in the introduction, we could 
in principle deconvolve out reconstruction with the known 
PSF. This, unfortunately, would completely remove the ef- 
fect of the Wiener filter leaving us with the extremely noisy 
and radially oscillating reconstruction of the no-prior esti- 
mator. 



4.4 Signal-to-noise of matter haloes 

Here we would like to compare the expected pixelised signal 
in a reconstruction to the expected noise level within the 
same pixel. 

If the density modes, 5(1), of Sect. 14.31 belong to a ran- 
dom field - such as the noise in the reconstruction or the 
random density fluctuations in a particular fiducial cosmo- 
logical model - the variance of the density contrast inside 
the a pixel will turn out to be 



3-D mapping with 3-D weak lensing 13 




0.1 0.20. 30. 40. 50. 60. 70. 80. 9 1 1.11.21.3 0.10. 20. 30. 40. 50. 60. 70. 80. 9 1 1.11.21.3 

redshift redshift 



Figure 7. The radial PSF of a fiducial survey as in Fig. [3] (ft = 30arcmin — 2 ; curves continue up to redshift z = 1.5). A density pixel 
of radius 1' sitting on a lens-plane with certain redshift (numbers attached to curves) is in the process of a reconstruction attenuated 
and spread out. The left figure assumes a transverse Wiener filter with a = 10 — 2 , whereas the right figure uses a radial Wiener filter 
with /3 = 10 — 3 . For example, a pixel at z = 0.05 is shifted to z ~ 0.2 (maximum of PSF) and attenuated with a factor of ~ 0.03 with 
transverse filtering (thick solid line). Note that the PSF of the radial filter can become negative. 
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Figure 8. These plots quantify the shape of the radial PSF for circular pixels of radius 1': the redshift of the maximum response (left 
panel) and the z-width (lc) of the PSF (right panel) as function of redshift of the original pixel. The shape of the PSF also depends on 
the tuning parameter a (different line styles). See Fig. [7J for a concrete example of the PSFs. The plots are based on a fiducial ACDM 
model with WMAP3-like parameters, an intrinsic source shape noise of <r e = 0.3, a source number density of n = 30arcmin~ 2 and a 
source redshift distribution as in Eg. 1701 The sources are binned between z = 0... 1.5 in equally spaced bins with Az = 0.05. To obtain 

the corresponding plots for different a' f or n', but all other parameters kept, one has to rescale q m q x M) \Wr see Sect. 14.11 



r(e 8 



i 

2~7 



Al£\F(£Q s )\ 2 P(£) 



(83) 



for a power spectrum, P(£), quantifying the fluctuations in 
5(1). This follows from Eq.[75l 

If we apply this integral transformation to the power 
spectra (matrices) Eqs. [62] and 1641 we can transform the 
Fourier space representation to signal-to-noise estimates in- 
side pixels as function of lens-plane redshift. From that one 
can infer that the average signal-to-noise inside a circular 
pixel will be approximately the signal-to-noise of one effec- 
tive £-mode, Fig. [3] As the integration kernel, x\F(x)\ 2 , in 
Eq. [83] peaks at about £Q B « 1.36, the effective mode fre- 
quency will be near £ e g « 4666 (0 s /l') _1 - 



The expected covariance of noise inside lens-plane pixels 
arranged along the line-of-sight will be: 



2tt 



Cov(5(e s ,a)) = — / d££\F(m s )\ 2 N M v(a) 



(84) 



This can be directly compared to the blurred signal from a 
single halo (Eq. [77]), sitting at some redshift, being filtered 
by our Wiener filter: 

g(&(3.,q))= , M& _ B ' a) ■ (85) 

Note that we only consider the smoothed signal about the 
halo centre here, where we expect to see the highest signal- 
to-noise of the halo reconstruction. The Wiener filtering does 
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Figure 9. Same as Fig. [8] but using a radial Wiener filter. For small (3 the curves disperse in the right panel because the radial PSF 
starts to oscillate about zero similar to the filter without prior. 



not bias the angular position of a halo (the phases of the 
Fourier modes), in particular the centre position. Therefore, 
we will, on average, find the halo centre in the reconstruction 
in the original direction as before filtering. 

For a singular isothermal sphere (SIS; e.g. 



(86) 



iBartelmann et aLll200ll ) we have 



5(6) 



47rc 2 a(ui h ) 



/k(wh) 



where ov is the velocity dispersion inside the SIS and a(wh) 
the scale factor at the radial comoving distance of the halo 
centre, w^. The Fourier transform of the SIS profile is 



S(£) = 



8ir 2 c 2 a(wh) 



(87) 



Aw V C / fK(w\C)£ 

As the signal inside a pixel, Eq. 1771 is linear in the 
original underlying amplitude of 5, whereas the noise inside 
the pixel is independent of S (Eq. I84[l we can already infer 
that the signal-to-noise of the SIS in the reconstruction has 
to scale with S/N oc a 2 , roughly the mass of the halo. We 
can expect this to be approximately true for more realistic 
halo profiles as well. 

More specifically, IBartelmann et all ||200jjT ) give for the 
virial massj of the SIS 



M 2 



where 



WGH(z) 

6.58 x 1O 14 M h' 1 ( 



10 3 kms 



E(z) 



(89) 



e(z) = y/n m (i + z) A + (i - n m - n A )(i + z) 2 + Ov (90) 

is the normalised Hubble parameter, H(z)/Hq, as function 
of redshift and G Newton's constant. For f2 m = 0.23 and 
z < 1, E(z) w 1 — 0.345 z provides a fair approximation. 
Therefore, the scaling for a SIS is S/N oc M^qq. 

Fig. [10] displays the signal-to-noise of the central pixel 



5 The mass contained in a radius at which the mean density of 
the SIS equals 200 times the critical density of the Universe. 



of a smoothed, reconstructed SIS. The calculations are done 
for a SIS with fixed velocity dispersion but can easily be 
rescaled for other a v , or different fiducial surveys with varied 
source density or shape noise. See the caption of the figure 
for details. 

The signal-to-noise of the central pixel is highest for 
small halo redshifts and decreases rapidly towards higher 
redshifts. This is because the number of sources behind the 
halo decreases as we move towards higher halo redshifts. 
Therefore, the technique identifies haloes most efficiently at 
low redshifts and becomes increasingly ineffective towards 
higher redshifts. We also found by varying the redshift dis- 
tribution of the sources, while keeping their number den- 
sity and intrinsic ellipticity distribution constant, that the 
signal-to-noise fall-off is shallower for a deeper survey, al- 
though the z-spread increases owing to the same number of 
sources being distributed over a larger z-range. 

We also gather from the results, that the signal-to-noise 
decreases as we tune to smaller values of a, thereby reduc- 
ing the Wiener regularisation. For a transverse filter and 
the given fiducial survey, values of a < 10 -2 are necessary 
to reduce the bias of the Wiener filter (lcr-z-spread and z- 
shift less than ~ 0.2) to an acceptable level, Fig. [S] For our 
radial filter, a value of f3 < 10 -3 should be used. We ex- 
cluded values smaller than a ^ 10 -3 (radial: f3 ^ 10 -4 ) 
as this brings us too close to the no-prior filter which ex- 
hibits heavy, undesired oscillations in radial direction and 
very poor signal-to-noise. Note that a also depends on the 
redshift distribution of the sources. 

An example reconstruction (transverse filter) of a set of 
SISs planted into a one-square-degree field-of-view between 
redshift z = 0...1 can be found in Fig. 1111 This example 
underlines the theoretical expectation for our fiducial sur- 
vey that the signal-to-noise of the haloes in the reconstruc- 
tion drops quickly if we move towards higher redshifts; the 
most distant haloes (F-I) are barely or not distinguishable 
from the background noise. Note that the signal-to-noise 
levels cited in the figure caption were computed based on 
500 FFT-reconstructed noise realisations, see Appendix lA3l 
One can also observe the effect of the PSF that spreads out 
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Figure 10. Signal-to-noise of a reconstructed SIS, positioned at different redshifts, x-axis, with a v = 10 3 kms — 1 ; this corresponds to 
-M200 = 6.6 X 10 14 Mq/i _1 (1.6 X 10 15 Mq/i -1 ) at 2 = (1). Plotted is as function of tuning parameter (line styles) the signal-to-noise at 
maximum response. The assumed fiducial survey has the parameters outlined in Fig. \E\ To obtain the corresponding plots for different 

a' e , a' v or n', but all other parameters kept, one has to rescale a 1— » a X f~"J ( Hi 1 an d § ^ H x (ct^) (o'y («') ,as ^^cribed 
in Sect. 14.11 Left figure: transverse Wiener filter. Right figure: radial Wiener filter. 




Figure 11. A set of nine SISs (<r v = 10 3 kms -1 ) with redshifts 2 = 0.1, 0.2, . . . , 0.9 (A-I in alphabetical order) has been planted into 
a particular noise realisation of the fiducial survey defined in Fig. [8] The left panel is a 2-D projection of the reconstruction onto the 
sky while the right panel is a projection showing the extension of the volume in redshift. The differently coloured contours are constant 
signal-to-noise levels corresponding to S/N = 8, 6, 5, 3, 2.5 (darkest to brightest colour). The reconstruction is smoothed with a kernel of 
radius 1'. The employed transverse Wiener filter has a = 0.01. Statistical errors can dissolve peaks as can be seen for D and E or shift 
peaks towards wrong redshifts, see D for instance. Note that all redshifts are biased to some extend, left panel of Fig. [8] 
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Figure 12. Reconstruction of the lensing convergence on a 
60' X 60' patch falsely assuming that all shear signal originates 
from matter inside the patch (left: B-mode, right: E-mode). Ac- 
tually, the signal is solely produced by a SIS, <r v = 10 3 kms — 1 
and z = 0.2, sitting in a separation of 6' off the right edge of the 
patch. The noise in the images is due to the finite number density 
of (intrinsically circular) sources, n = 30arcmin — 2 , used for the 
fiducial survey here. 



the peaked matter distributions and moves, for instance, 
the maximum of the response of low-redshift haloes towards 
higher redshifts. Due to the noise in the data, some peaks (C 
and D) fragment into different parts and turn up at clearly 
wrong redshifts. An analogous reconstruction (/3 ~ 10 -3 ) 
employing a radial filter looks essentially the same. 



4.5 Edge effects 

The reconstruction algorithm assumes that the cosmic shear 
signal within the patch of the survey originates solely 
from matter within the field-of-view. Thi s is the same for 
all K aiser-Squires type algorithms (KS; iKaiser fc Squires! 
1 19931 ). This introduces, especially at the boundaries of the 
field, a bias if the actual source of the shear within the field 
is outside the patch. This problem would not exist for a 
full-sky implementation of the algorithm. 

To have an qualitative assessment of the impact of this 
bias, we consider one particular, extreme case which is shown 
in Fig. Q2] For this mock survey, a square patch with 60' on 
each side and the previous redshift distribution of sources, 
we made up the case that there are no sources of shear at all 
inside the patch area. The total shear signal observed stems 
from one massive singular isothermal sphere located close to 
the right edge of the patch. The figure shows a KS recon- 
struction (left: B-mode, right: E-mode) of the convergence 
inside the patch ignoring an outside source as possible origin 
of the shear. The 3-D reconstruction of the density contrast 
for a lens plane at about the redshift of the SIS exhibits 
exactly the same pattern. Lens planes well separated from 
the SIS redshift show no or little reconstructed mass den- 
sity while neighboured lens planes may be affected due to 
the radial spread of the reconstruction PSF. 

It can be seen from this test that a source outside the 
observation area can indeed produce mass density inside the 
reconstruction patch. However, this "ghost signal" seems to 
be mainly focused on the edges. The largest inferences can 
be seen on the edge closest to the SIS, namely at the middle 



and at the corners, and on the edge opposite to that with 
roughly the same but inverted signal pattern. 

Furthermore, it can be seen that this ghost signal is 
associated with a B-mode signal of roughly the same am- 
plitude than the E-mode signal, also most prominent at the 
edges (closest and opposite) but leaking somewhat more into 
the patch. Therefore, if an off-patch source produces a sig- 
nificant signal inside the patch area it may reveal itself by 
an equally significant B-mode signal. 

Note that the response pattern looks different if we use 
FFTs for the reconstruction, as they assume periodic bound- 
ary conditions. With FFT the SIS can actually be recon- 
structed as mirror image near the edge opposite to the edge 
closest to the true position of the SIS. Yet, we find that even 
with periodic boundary conditions a B-mode signal is gener- 
ated inside the patch which has the same order of magnitude 
than the E-mode signal. 

This bias can be controlled by making the reconstruc- 
tion area larger than the patch area, attributing to this addi- 
tional patch frame infinite noise. By doing this, we explicitly 
allow the reconstruction algorithm to place sources of shear 
outside the observed patch area. 



4.6 Reduced shear 

Throughout this paper we assumed the weak lensing ap- 
proximation to be perfectly valid, i.e. ellipticities of galaxy 
images are unbiased estimators of the cosmic shear and not 
the reduced shear, y/(l — ft), as they actually are. 

If the distinction between reduced shear and shear is 
ignored, the convergence will in the vicinity of density peaks 
biased towards smaller values and hence the density will be 
biased towards smaller values as well. 

Presumably, the algorithm can be modified in this re- 
spect by running the linear algorithm on (1 — Kj ), where 
is k = QS is the minimum variance convergence reconstruc- 
tion of the previous run; for the first run one starts off with 

= 0. For a 2-D convergence reconstruction this iterative 
approach is k nown to converge quickly to the right solution 
(Scliak 1998, Sect. 4.3). The noise of the nonlinear recon- 
struction would have to be estimated by a series of noise 
realisations. We will postpone a test of this approach for a 
3-D reconstruction to a future paper. 



5 OUTLOOK 

For a demonstration and an outlook of what may be achieved 
with this technique in the near future we apply the recon- 
struction method to a simulated lensing survey patch of four 
square degree size. The simulated survey has a source red- 
shift distribution as in Eq. [70] but now with zo = 1.0, a 
source number density of n = lOOarcmin -2 , corresponding 
to 1.6 x 10 6 sources, and an average intrinsic ellipticity vari- 
ance of <7 e = 0.3 for all sources. It is further assumed that we 
have (photometric) redshifts for all individual galaxies and 
that we are unable to obtain (photometric) redshifts beyond 
a redshift of z = 2, whence we truncate p(z) at a redshift of 
two. The mean redshift of the sources is therefore z = 1.1. 
These survey parameters roughly reflect a space-based lens- 
ing survey which we consider as the best possible choice for 
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Figure 13. 3-D mass density reconstruction of a 2deg X 2deg patch of a mock lensing survey with lOOarcmin -2 source number density 
and mean source redshift z = 1.1 with z $C 2.0. The patch was chosen to have a large number of high-density peaks. The intrinsic shape 
variance is chosen to be <j e = 0.3. The reconstruction was done on a 128 X 128 pixel 2 grid with 25 source z-bins (Az = 0.08), and 20 lens 
planes (Az = 0.1). The reconstruction was done on grid larger than the observed area, allowing for a frame with 7'.5 width. A transverse 
Wiener filter with a = 0.05 was employed. For the signal-to-noise maps, which are shown here, the data was randomised 100 times. 
The upper left (projection on the sky; x vs y), lower left (x vs. redshift) and lower right panel (y vs. redshift) show the signal-to-noise 
data cube of the reconstruction with 3-D iso-levels of S/N = 6, 5, 4.5, 4, 3.5, 3, 2.5. The upper right panel shows for comparison true high 
over-densities (not over-densities due to projection) on the lens-planes used for the ray-tracing. The most significant peaks, S/N 4, 
can be identified and matched with true matter peaks (solid circles). Several less significant peaks can be matched as well (not all 
encircled). But there are also peaks in the reconstruction with no apparent real counterpart (dotted circles) or high density peaks in the 
true distribution with no counterpart in the reconstruction (dashed circles); not all cases have been highlighted. The numbers denote the 
estimated (left) or true redshifts (right) of the structures. The most significant peaks in the upper left panel are labelled with alphabetic 
letters, A-G, and can be identified in the radial projections in the lower row. Label X denotes a case discussed in the text. No attempt 
has been made to correct for the z-shift bias which is most evident for z < 0.15. 



mass reconstructions due to the small ratio a e /\fn, which is 
the main factor in the signal-to-noise of the maps. 

The mock survey was generated by ray-tracing through 
an TV-body simulati on run with t he publicly available ver- 
sion of GADGET- 2 (|Springel2005l ). We have used 256 3 Dark 
Matter particles in a box with a side length of 150 Mpc, 



which together with the adopted ACDM cosmology (fi m = 
0.25, fi b = 0.04, Q,a = 0.75, a g = 0.78) results in a particle 
mass of 1.2 x 10 10 h' 1 M©. We have traced 2048 2 light rays 
trough 25 matter slices with a thickness of 150 h Mpc up 
to z = 2 using the stan dard multiple-lens-plane algorithm 
(e.g. Irlilbert et al]|2009l . and references therein). This yields 
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the Jacobian matrix of the lens mapping to the back side 
of each matter slice. We have then created the mock survey 
by randomly sampling the slice boundaries with galaxies ac- 
cording to the desired redshift distribution, onto which we 
interpolate the shear from the light ray positions. 

For the reconstruction, we subdivide, equally in red- 
shift, the source catalogue into 25 sub-samples with width 
Az = 0.08. This number of source bins is fixed here by 
the number of lens planes originally used for ray-tracing to 
obtain the simulated data. In reality, Az reflects the aver- 
age uncertainty in the redshift estimates which is probably 
somewhat too pessimistic with 0.08. However, we would like 
to stress that a much finer binning does not make much 
difference in the reconstruction - at least if the redshift esti- 
mates are not biased and the redshift distribution inside the 
bins is accurately known - since the lensing efficiency is but 
a slowly changing function with source redshift. This is re- 
lated to the known fact that increasing the number of source 
bins beyond a few, in constraining cosmological parameters 
with lensing tomography, does not significantly increase the 
constraints. Hence, we do not expect a notable improvement 
by using a larger number of source bins. 

The number of lens planes, on which the matter density 
contrast is to be estimated, is set to twenty ranging between 
z = . . . 2 with Az — 0.1. We do not use more lens planes, 
as a radial resolution exceeding Az — 0.1 should not be ex- 
pected according to the foregoing discussion. In the following 
plots, we interpolated between the lens planes. 

Fig. 1 13 1 shows the reconstruction of the simulated patch 
as signal-to-noise map either projected onto the sky or in two 
different radial projections. Note that the sky-projections 
are actually 3-D signal-to-noise maps. A signal-to-noise map 
of the projected matter distribution, as commonly extracted 
from lensing surveys (smoothed convergence maps), can be 
anticipated to have higher significance. For comparison also 
a sky-projection of the most massive density peaks on the 
original lens planes is displayed. For this demonstration a 
patch particularly rich in over-density peaks was chosen, on 
average the number of peaks in a four square degree patch 
in a WMAP3-like universe should be expected to be smaller. 

We observe that although we have a quite optimistic 
lensing survey, the significant features in the reconstruction 
mainly correspond to the most massive haloes in the Dark 
Matter density field. A detection, S/N > 3, is restricted 
to the regime z < 1. Below S/N w 3.5 we find cases in the 
reconstruction which do not seem to have a real counterpart. 
Those can be explained as pure statistical flukes of the noise 
pattern which are still to be expected on this signal-to-noise 
level. It may be that 100 noise realisations, as performed 
here, do not suffice to estimate the noise pattern accurately 
enough everywhere. We also seem to have cases where we 
have no significant detection of arguably high peaks in the 
true Dark Matter density. The problem may be here that 
we are locally lacking the required sampling of background 
sources. Peaks along the same or close to the same line- 
of-sight are not always distinguished from each other and 
merged into a mediate redshift instead (see E), although 
there are cases where the algorithm succeeds in doing so 
(see C). 

Quite interesting is the case X, upper right quadrant, 
where we have a large peak at very low redshift z = 0.02 
which however is completely missed by the reconstruction or 



might be hinted at by a very weak detection at a completely 
wrong redshift. The explanation for this case probably is 
that the lensing efficiency for a structure at a redshift that 
low is either too small for sources at higher z or the num- 
ber density of sources which in principle would be sensitive 
enough to light-deflections at that low redshifts is too small 
(low z). In this sense, the lensing approach appears to be 
blind towards structure at very low z. 

Seemingly, the fiducial survey does not perform a lot 
better than the inferior survey in Fig. 1111 where we assumed 
n = 30arcmin -2 and a somewhat shallower mean redshift. 
An increase in n from the latter to the former survey, al- 
though observationally rather challenging, only improves the 
average signal-to-noise by a mere factor of ^1073 ~ 1.8. 
This again highlights the problem of all lensing mapping 
schemes. 

The signal-to-noise of the reconstruction can be 
strongly increased by using a more conservative, larger, tun- 
ing parameter (here: a — 0.05). For instance, structure D 
can be boosted to a signal-to-noise of approximately ten by 
using a = 1, however on the expense of loosing information 
on the redshifts, revealed by a significantly larger spread in 
radial direction, and more z-shift bias. 

Note that the z-shift bias is somewhat different in a 
signal-to-noise data cube as opposed to a density contrast 
data cube which was discussed in the foregoing section. 
Moreover, we have not tried to correct this bias in our red- 
shift estimates. In the plot, therefore, particularly structures 
at low redshift, z < 0.15, are artificially shifted to larger red- 
shifts, z ~ 0.2. 

The redshift estimates of structures more significant 
than S/N ~ 4 are all in all reliable within ±0.1, provided 
they are not unresolved structures at different redshifts, but 
become increasingly inaccurate below that signal-to-noise 
level. 

The scenario discussed here contemplates a case where 
we have a small but deep patch of which we would like to 
obtain a 3-D map with a preferably good spatial resolu- 
tion. Another scenario would be a (almost) full-sky survey 
to which heavy smoothing is applied to study the very large- 
scale density modes. Fig. [3] implies that for those modes the 
reconstruction may become easier than for the modes on a 
small patch as the signal-to-noise improves for lower I and 
the z-shift bias becomes less severe. 

The situation may be further improved if the shear 
data is combined with higher-order lensing information, such 
as gr avitational flexion, that may be avail able in the fu- 
ture dBacon et all 120061 ; lOkura et ail 120071 ). In the weak 
lensing regime, higher-order moments of lensing distortions 
can be modelled as linear functions of the lensing con- 
vergence, as derivatives of «. Therefore, higher-order mo- 
ments are easily implemented into the algorithm as new ex- 
tra set of redshift bins with a new linear projection opera- 
tor "Pf k " as opposed to P 7K , Q remains unchanged. Very 
crudely - same number of sources as for cosmic shear and 
same (uncorrelated) signal-to-noise provided - this will be 
equivalent to effectively increasing the number density of 
sources in a conventional cosmic shear survey by a factor 
of N + 1, where N is the number of new lensing distortion 
moments that is included; the signal-to-noise in the recon- 
struction consequently goes up by ViV + 1. Actually, as al- 
ready known from first 2-D applications (|Bacon et al.ll2006l ; 
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iLeonard et al.l [20071, the 3-D maps will become better on 
small angular scales with less improvement on larger scales 
as the signal-to-noise of flexion- sampled mass density fields 
is scale-dependent with most information on small scales. 
The radial resolution of the maps, however, will only mod- 
erately benefit from the new information, i.e. no more than 
from an extra cosmic shear catalogue with same size, as the 
dependence of higher-order lensing on redshift is identical to 
that of cosmic shear. 



6 SUMMARY AND CONCLUSIONS 

In this paper we presented a linear algorithm - and two 
practical ways of implementing it (Appendix I A 1 1 A2|l - that 
allows to reconstruct the three dimensional distribution of 
matter based on a lensing survey with redshift information. 
All available redshift information, accurate source redshifts 
binned into thin slices and broad redshift information such 
as wide redshift distributions for faint, high-z source sub- 
samples, can be combined to find a matter distribution that 
should be on average closest to the true distribution. The 
statistical uncertainties of the source redshifts, attached as 
PDF's of the redshifts to the corresponding source subsam- 
ple, and the statistical uncertainty of the shape measure- 
ment can be properly factored into the reconstruction. Prin- 
cipally, as shown in the paper, it is also possible to account 
for intrinsic alignments of the sources, that have an impact 
on the expected shape noise covariance entering the filter, 
and shear-intrinsic alignments. 

The presented alg orithm, being a gene ralisation of the 
algorithm proposed bv lHu fc Keetonl (|2002h . yields the most 
probable matter distribution, if the matter density field 
obeys Gaussian statistics, or a minimum-variance solution 
otherwise. Due to the heavy smoothing required in a 3-D 
matter distribution reconstruction based on lensing tomog- 
raphy, the reconstruction is probably close to the maximum- 
likelihood solution. 

The algorithm is a Wiener filtering of the lensing to- 
mography data. The strength of the Wiener filter prior, i.e. 
the assumed second-order correlations between densities in 
the reconstruction, can be tuned with an additional param- 
eter, a. We looked into two different cases for a regularisa- 
tion, namely one case for which only correlations between 
densities at the same redshift are given (transverse filtering) 
and a case where only correlations along the line-of-sight are 
specified (radial filter). Both priors show similar properties, 
it is not possible to single out one filter as superior filter. 

We showed that a reconstruction with virtually no prior, 
a ~ 0, gives an unbiased albeit extremely noisy reconstruc- 
tion. Moreover, a reconstruction with small a exhibits strong 
radial oscillations. The oscillations vanish and the signal-to- 
noise improves considerably when one increases a, thereby 
enhancing the effect of the Wiener prior. However, a Wiener 
filtered reconstruction represents, as all Wiener filter esti- 
mators, a biased image of the true underlying signal. This 
bias can in our context be expressed as a radial PSF or re- 
sponse of the reconstruction for which we give an analytic 
expression in Eq. [HTJ The effect of the PSF is that a peak 
in the original density field will, on average, be stretched 
out and shifted in radial direction, see Fig. [7] The stretch 
and shift depend on the fiducial parameters of the survey 



and the tuning of the filter. The bias essentially reflects our 
inability to pin down exactly the redshift of a density peak 
due to the noise in the data. 

Generally, the larger the a-tuning, the better the signal- 
to-noise of the reconstruction but the more biased be- 
comes the reconstruction. In the extreme case of over- 
regularisation, one obtains essentially a 2-D reconstruction 
on the sky stretched out over the entire radial range lacking 
any radial information, providing the best signal-to-noise of 
a map of the projected matter distribution though. 

As strategy, we suggest to choose a parameter a with 
moderate z-shift bias, see e.g. Figs. [8] and [9] and to relabel 
the redshifts of the reconstruction lens-planes according to 
the expected statistical average of the z-shift bias. At low 
redshift the radial PSF flattens out, as can be seen in these 
figures (there: z < 0.3). A similar behaviour can be observed 
at the high redshift end of the survey. In these regions, we are 
hence unable to uniquely relabel the lens plane redshifts. At 
the low redshift range of that fiducial survey, all we can say 
is that the structure generating a peak in the reconstruction 
has to be located somewhere below z < 0.3. 

The reconstruction will be relatively noisy so that a fi- 
nal smoothing in 2-D (lens planes only) or 3-D, for example 
with a Gaussian kernel, has to be applied. The thereby ob- 
tained smoothed reconstruction is still a minimum-variance 
reconstruction of the true density field subject to the same 
smoothing, an invertible smoothing provided. 

As can be seen from the implementation strategy, the 
algorithm is in essence a series of linear operations - pro- 
jections, pixel-rescalings and convolutions - that have to be 
applied step by step to the pixelised input data. We would 
like to point out here that for the derivation of this algorithm 
we can equally imagine the data to be gridded on spherical 
shells in a spherical coordinate system; Fourier modes would 
be spherical harmonics on the unity sphere instead of waves 
on a plane. In fact, by doing so we will get exactly the same 
algorithm that, however, now tells us to linearly combine 
different lens-shells or to perform convolutions on the unity 
sphere. Therefore, the recipe described in Sect. lAll is a recipe 
for a full-sky reconstruction as well, for which the discussed 
optimisations and biases apply too. 

We also gave an analytic estimate for the noise covari- 
ance of density pixels in the reconstruction, Eq. 1841 for an 
idealised survey with homogeneous noise, no gaps and infi- 
nite extent. Several numerical algorithms for the covariance 
of a realistic survey are outlined in Appendix IA3I A com- 
parison of the noise to the signal of a SIS, Eq. 1771 shows 
that the signal-to-noise of a matter halo quickly decreases 
towards increasing redshifts in a reconstruction (Fig. I10[) . 

We noted that the algorithm (if implemented on a flat 
sky), as other similar reconstruction methods known in lens- 
ing, implicitly assumes that all cosmic shear signal originates 
from matter fluctuations inside the reconstruction area. If a 
shear of galaxy images inside the patch actually stems from 
a source outside the patch, this will introduce ghost images 
near the edges of the patch. In principle, the generation of 
significant ghost images can be detected by significant fea- 
tures in the mass density reconstruction based on the B- 
modes of the shear field, which come along with the E-mode 
ghost structures, and have the same amplitudes. These edge 
effects can be reduced by putting an additional frame around 
the reconstruction patch with infinite noise. 
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A fiducial survey with 30 sources per arcmin 2 , mean 
redshift of z = 0.85, a e = 0.3 is unable (S/N < 3) to identify 
a matter halo with mass ~ 7 x 10 14 Ii~ 1 Mq or less beyond 
redshift z ~ 0.6 in a reconstruction with moderate (a ~ 
10~ 2 , transverse filter) z-shift bias. At a lower redshift of 
z ~ 0.1 the detection limit corresponds to a mass of ~ 1 x 
10 14 fickle. 

This redshift limit increases for deeper surveys which 
we demonstrated by doing a reconstruction of a mock sur- 
vey that broadly mimics the characteristics of a future space- 
based lensing survey. Still, even boosting the number density 
of sources to 100 arcmin -2 and going deeper to z = 1.1 only 
moderately improves the performance of the reconstruction. 
For such a survey, structures in the map of S/N w 4 or 
greater have accurate redshifts within ±0.1, below this level 
we have to expect false positives. Generally, we have to ex- 
pect to miss out on higher density peaks, even for a powerful 
lensing survey as the one assumed here. 

The situation can be expected to become better if a full- 
sky reconstruction with heavy smoothing is done to study 
- in contrast to structures on non-linear scales - density 
modes on very large scales: the signal-to-noise increases to- 
wards smaller £, the bias effects in the reconstruction ease. 
We expect the inclusion of higher-order lensing distortions 
into the 3-D reconstruction to further improve the recon- 
struction, very roughly increasing the effective number of 
"cosmic shear sources" by a factor of N + 1 for N extra 
distortion moments. This, however, needs to be explored in 
more detail elsewhere. 

The authors see cosmography as main application of the 
reconstruction technique. With 3-D lensing conventional 2-D 
lensing convergence maps can be supplemented to some de- 
gree by redshift information, mainly restricted to larger mass 
overdensities or to large density modes. This may not give us 
the full 3-D appearance of the mass distribution of, say, indi- 
vidual galaxy clusters - rather the 2-D projection stretched 
out in radial direction about the maximum-likelihood ra- 
dial distance - but may enable us to disentangle physically 
unconnected structures along the line-of-sight providing a 
better understanding of the matter distribution inside a 3-D 
volume. This can be vital for a more sophisticated modelling 
later on. 

As an obvious application, the 3-D data cubes may be 
used to make catalogues of matter overdensities endowed 
with redshifts without any assumptions about density pro- 
files. Estimating the density profiles of individual structures, 
essentially projected to 2-D, in the 3-D lensing map is con- 
ceivable, too, but only if the full reconstruction PSF, includ- 
ing the transverse direction, is accounted for. Owing to the 
map biases an estimate for the mass of structures should 
be obtained by directly fitting a density model to the shear 
data once they have been identified in the non-parametric 
3-D reconstruction. 

In principle the 3-D maps can be compared to the 3- 
D distribution of galaxies as derived from galaxy redshift 
surveys. In this case, again, care has to be taken for the bi- 
ases and limited radial resolution in the lensing map. For a 
fair comparison, this may be done by convolving the map 
of galaxies with the radial PSF of the lensing maps. For a 
statistical comparison of galaxy and (lensing) matter maps 
requi red for studying the galaxy biasing (e.g. ISimon et al.l 
120071 ). however, we would suggest to forward and fit a 3- 



D model of matter density fluctuations and galaxy biasing 
to the 3-D shear tomography and its cross-correlation with 
galaxy positions, rather than producing a biased Wiener fil- 
tered 3-D matter map that is cross-correlated with the 3-D 
distribution of galaxies. In the best case, with all biases in 
the 3-D mass map being fully understood, both ways should 
give same results anyway. 

The same holds for attempts to estimate matter density 
power spectra from the 3-D lensing maps. As shown in this 
paper, due to low signal-to-noise the power spectra on the 
lens planes will be strongly influenced by the Wiener prior 
which has to be taken into consideration for an estimate 
of the density fluctuation power. Even if this is done prop- 
erly, it should only give a result equivalent to fitting a 3-D 
model for the shear tomography c orrelation functions which 
is straightforward in comparison (Mas sev et alj |2007) . 
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APPENDIX A: OPTIMISED 
IMPLEMENTATIONS OF THE ALGORITHM 

Al Real space reconstruction 

As seen in Sect. 13.21 finding the MV-solution for some given 
data vector of shear tomography ellipticities requires, in the 
case of the density contrast, the evaluation of 

<5mv = [la + SiQ^N^P^Qn'SiQ'P^N-S . (Al) 

T his filter may be used for a reconstruction. As pointed out 
bv lTegmarkl (|l997T ). the estimator Eq. [7]is identical (a / 0) 
to 

s M v = SR^RSR 1 +aU]- 1 d , (A2) 

which is a form that is also often encountered in the liter- 
ature. For that reason, we can equally write as linear filter 
for the density contrast field as: 

Smv = (A3) 
S f Q*P^ [NJ^QSgQ*P^ + aN^No + al]" 1 IM^S . 

stc P 3 =M - \«rtep 2 stcp 1 

Casting the estimator into this form is actually more conve- 
nient from a numerical point of view if we have to consider a 
non-diagonal noise covariance, N, e.g. in the case of intrinsic 
alignments. Here, we split the noise covariance, N = Nd + N , 
into the sum of diagonal, Nd, and off-diagonal, N , elements; 
the diagonal matrix Nd is easily inverted and multiplied with 
a vector. 

For an implementation of this estimator we proceed 
in three steps. For the first step, see the underbraces, the 
binned shear data 7 is multiplied by the inverse of the diag- 
onal elements of N. This yields the intermediate result "a;" 
that is actually a set of shear grids in which each grid point 
is weighed by its noise variance. For the second step, x has 
to be multiplied by the matrix M 1 . We employ the method 
of conjugate gradients to invert the equation My = x with 
respect to y (|Press et alJll992h . This essentially boils the 
matrix inversion down to a series of multiplications with M 
until the solution converges. Again, multiplying by the ma- 
trix M is treated as a one-by-one multiplication by the series 
of individual matrices contained within the square brackets. 
As initial guess for y for the iterative algorithm one com- 
monly uses y — which appears to work fine in the context 
of this paper. Finally, for the last step y is multiplied by 
SsQ'P^. 

Multiplying a matrix with a vector means the applica- 
tion of a linear operation on the input vector. Here, each 
matrix has a different effect on the vector that is hard-wired 
as vector- in/vector-out sub-routine individually. The multi- 
plication with Q (or its transpose Q') is easily implemented 
as it is just summing J]\ QijX k along every line-of-sight Ok- 
Multiplications with inverses of diagonal matrices are also 
easy as mentioned above. 

The time consuming matrix operations are connected 
to P 7 k, S«5 and N , if considered, which denote convolu- 
tions of the input vector. Only in the situation where we 
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are using radial Wiener filtering, S$x will become an easy 
operation requiring only sums along line-of-sights analogous 
to Qx. The most effective way of doing convolutions uses 
Fast-Fourier- Transforms (F FTs) as an efficien t application 
of the convolution theorem |Press et al.lll992T ). For that we 
need to FFT the different layers of x (either lens planes or 
source bins) separately and multiply the Fourier coefficients 
with the power spectrum of Ss, the Kaiser-Squires kernel 
D(l), in case of P 7K , or the power of intrinsic alignments; 
the power spectra are functions of lens planes or source bins. 
The window function of the grid cells has to be taken into 
account as well. The transpose of P 7 „ corresponds to D* (£) . 
After that we re-FFT to real space. Note that the whole 
chain of operations P^QSaQ*?!^ can be done completely 
in Fourier space with no back and forth FFTing in between. 

The FFT method has two drawbacks: i) It requires grids 
sizes which are powers of two so that not any arbitrary field- 
geometry can be processed, ii) Also, as well known, we will 
have edge effects as FFTs are assuming periodic boundary 
conditions (aliasing) . Regarding i) we have to find the small- 
est bounding box to enclose the whole field of view at the 
desired resolution. Resolutionwise a number of grid cells of 
128 or 64 for each coordinate axis is usually enough. With 
respect to ii) we can use the zero-padding technique or down- 
weighting towards the edges (see next section) to mitigate 
the aliasing. 

There is an alternative technique known in the liter- 
ature that speeds up convolutio ns but works solely in real 
space (|Padmanabhan et al.ll2003h , does not suffer from FFT- 
aliasing and the geometry constraint and, in fact, does not 
require griding. This approach is a real alternative to FFT 
but considerably slower as we found. As both techniques 
yield essentially the same reconstructions, apart from the 
edges, we propose to stick to the FFT technique in this con- 
text. 

Processing 25 lens planes and 40 source redshift bins, 
all with a grid size of 128 x 128 pixel 2 (the data vector, com- 
plex, and matter density vector, complex (E-, B-mode com- 
bined) , have sizes of 1.3 x 10 6 and 8.2 x 10 s elements, respec- 
tively) requires for the algorithm ~ 60 MByte of computer 
memory (double precession). The computation time for the 
second step on a AMD Athlon 64 processor with 2.4GHz 
is approximately five seconds for every iteration of step 2, 
which is the main bottleneck of the algorithm. For noisy 
data and (full, i.e. a = 1) Wiener filtering usually about five 
iteration steps are needed, whereas a reconstruction with 
no prior or heavily tuned-down Wiener filter (a ~ 0.01) 
requires up to ~ 150 iteration steps, depending on the de- 
sired accuracy, because less regularisation is applied to the 
data so that more details are fitted. The benchmark param- 
eters used here are usually, with contemporary real data, 
unnecessary. For that 64 2 -grids, Ni p ~ 10 lens planes and 
N z ~ 20 source bins are absolutely sufficient bringing down 
the computation time of one full reconstruction to roughly 
ten seconds. 

After the algorithm has converged, we run a final (trans- 
verse) smoothing of the lens-plane grids with a Gaussian 
kernel. The smoothing is usually required as a down-tuned 
Wiener filter is applied that produces noisy images. The 
blurring can be done efficiently by again employing the FFT 
convolution or a kdtree implementation. 



A2 Fourier space reconstruction 

If an idealised survey is a good assumption for particular 
data, the Fourier space filter Eq.[5l]is, due to the small ma- 
trices involved, a very fast and efficient way to perform a 
3-D reconstruction completely in Fourier space. In particu- 
lar, the matrix inversions can be done explicitly. To apply 
this filter, we FFT the gridded source z-bins, arrange the 
shear Fourier coefficients for each I inside one vector, 7, 
and perform Eq. [51] for every angular mode. The matrices 
involved have only sizes of the order of N z or/and Ni p , which 
are easily and quickly inverted. The Fourier reconstruction 
is eventually back-FFTed to acquire the real space mass den- 
sity field. Finally, the lens-planes need to be smoothed, since 
we are usually using a tuned down Wiener filter (a = 0.01, 
for instance). This smoothing can be done in real space or, 
more efficiently, directly in Fourier space immediately after 
the Fourier space reconstruction. 

Usually for I 7^ a regularisation (A 7^ 0), especially 
for the no-prior MV-estimator, is not necessary in a Fourier 
space reconstruction because the reason for the singularity of 
the filter is the mass-sheet degeneracy that affects solely the 
I = modes. Here, we lift this degeneracy by setting all I — 
Fourier coefficients to zero which amounts to saying that the 
mean density contrast on each lens plane is exactly zero. 
For a large enough field-of-view this should be a reasonable 
assumption. 

In order to reduce the FFT edge effects, the real 
space grids should, initially, be downscaled by sin (-7rA/2Ao) 
within a strip of width Ao, where A is the shortest distance 
of a pixel to the grid edge. 

The noise covariance is estimated by means of the well- 
known shot-noise formula [cr^ ] 2 jfii (Sect. [3]3}, ideally as- 
suming that the grid cells are homogeneously filled with 
sources of roughly equal intrinsic ellipticity distribution. 

A complete reconstruction on a Fourier grid for the 
same parameters (but 128 x 128 grids) as in the foregoing 
section takes less than a second. 



A3 Covariance of statistical uncertainties 

An estimation of the statistical errors, or more general the 
covariance, of pixel densities in the reconstruction is needed 
to convert the density reconstruction into a signal-to-noise 
3-D map. As a lowest-order estimate we suggest to use the 
analytical solution, Eq. 1841 that is based on homogeneous 
noise of the sources. It offers the r.m.s.-error for each lens- 
plane and the correlation of errors along each line-of-sight. 

As a final smoothing of each lens-plane will usually be 
necessary, the estimate has to be rectified slightly. Lets us 
assume that the errors, as, of the unsmoothed map are un- 
corrected, which will be roughly correct for a strongly down- 
tuned (small a) filter. By smoothing a lens-plane with a nor- 
malised kernel, K(x), one will combine an effective number 
of 



[J R2 d 2 0AT|fl|)f 

J R2 d^e[K(\e\)} 2 



(A4) 



pixels with uncorrelated errors as- The combined error of 
the weighted average (every smoothed pixel) will therefore 
be ag I viVeff . For a Gaussian kernel with width a one finds 
N c s — ty 2 a 2 . The truncation of the kernel due to the edges of 
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the map can be accounted for using the formula for N e g and 
setting K(x) = for areas outside the map. As the pixels 
one is smoothed over are slightly correlated, crs/V N c a will 
give only a lower limit for the true error (upper limit for the 
signal-to- noise), whereas as would be the noise upper limit 
being reached if all pixels were 100% correlated. 

In reality, one encounters surveys with varying source 
number densities and intrinsic source shot noise. One may 
desire to get a better estimate for the signal-to-noise of a 
pixel than provided by the foregoing formula, which is given 
for all pixels by 



wx 1 w f , 



(A5) 



where X" 1 = Q'P^N^P^Q, W = FtS«[X-% + al]" 1 
and F is the smoothing matrix (unity matrix in absence 
of final smoothing). However, the full covariance imposes a 
huge, if not almost impossible, computational task because 
of the enormous size of the involved matrices. 

What can easily be done, though, is to compute single 
matrix elements of the covariance: 



[N 



Vi 



v]N~ 
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(A6) 
(A7) 



where F, is defined as the ith column of the matrix F, which 
is the smoothing weight of any grid pixel relative to a fixed 
grid pixel i. This vector is filled with many zero because we 
are only smoothing within the same lens-plane. As indicated, 
the recipe is to compute the vectors Vi/j in a first step and 
then v\N~ 1 Vj in a second step. In order to attach error bars 
to single Smv one requires only one vector Vi = Vj. 

The definition of v t /j is almost identical to <5mv in Eq. 
lAll so that the algorithm discussed in Sect. I All applies with- 
out much change. On the other hand, this means the effort 
to obtain a matrix element of Nmv is roughly twice the ef- 
fort (Only off-diagonals; the same effort for diagonals) to 
remake a whole 3-D reconstruction, so we may only wish to 
work out the exact signal-to-noise or correlation of errors for 
a selected small number of pixels such as the peaks in the 
mass map that are suspected to belong to galaxy clusters. 

A third and most effective alternative for obtaining a 
signal-to-noise map is to randomise the ellipticities of the 
sources and to make a full reconstruction. With the algo- 
rithm outlined in Sect. IA1I this can be done hundreds of 
times in a reasonable time, maybe running parallelly on dif- 
ferent computers. The variance in the density contrast be- 
tween the different noise realisations would serve as estima- 
tor for the noise level in the reconstruction. 



how it can, in principle, be included by extra terms in the 
minimum- variance estimator for s. 

For that we assume that the signal/noise-correlations 
(here: correlations between the lens plane matter density 
and the intrinsic ellipticities of sources) between the various 
source z-bins are known either from measurement or from a 
physical model, U = (sn^). The minimum-variance Ansatz 
seeks to find a linear transformation, H, that minimises the 
residual matrix 



(Hd-s)(Hd-s) f 



(Bl) 



for all possible noise and signal configurations. For the sim- 
pler case that U vanishes, this residual matrix is, leaving out 
terms independent of H, on average 

SRtH 1 , 



H(RSR t + N)H t HRS 

which is minimised for 
H = SR^RSR 1 + Nl" 1 



(B2) 



(B3) 



If we now drop the condition of a vanishing U, the variance 
IB1I becomes more generally 

H(RSR t +N+RU+U t R t )H t H(RS+U f ) (SR t +U)H t (B4) 
If we substitute in the last expression 
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we recover exactly the variance IB2I with the already known 
minimising solution Eq. lB3l Undoing the former substitution 
in lB3l (the substitution is independent of H) gives the more 
general minimum- variance that now also accounts for signal- 
noise correlations: 

H = (SR 1 " + v^ll) [rSR 1 " + qN + a/qJRU + aAi/R 1 "] ' .(B7) 

Here, we already included the rescaling of Sw a _1 S 
needed for the Q-tuning (Saskatoon filter) supplemented by 
U i — > y/a 1 U to be consistent with the a-tuning (the signal- 
noise correlations are unchanged by changing a). We find 
that including the signal-noise correlation adds three new 
terms to Eq. IB3I In practice, this estimator can be tackled 
by the same numerical techniques as outlined in Appendix 

El 



This paper has been typeset from a TpjX/ BTjrjX file prepared 
by the author. 



APPENDIX B: SHEAR-INTRINSIC 
ALIGNMENT CORRELATIONS 



The estimator Eq. IA2I is the minimum variance solution, 
or the most likely solution for s in the Bayesian sense for 
pure Gaussian statistics, for a given data vector d pro- 
vided that there is no correlation between noise and sig- 
nal, i.e. (sn.^) = 0. In the context of gravitational lensing, 
however, these correlations can occur if, for example, the 
shear si gnal itself is correlated to the intrinsic alignment of 
sources i|Hirata fc Seliakll2004 ). We ignore this effect for 3-D 



mass reconstructions in this paper but would like to sketch 
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